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§ ■ ABSTRACT 

1—5 ; 

CN ! We propose a bipolar jet model for supernova polarization. Light from the 

,__! i main component of the supernova (which we call the bulk supernova for short) 

^ ! scatters off electrons in the jets and is polarized: this polarized light is added 

Tj- ! to the direct emission from the bulk supernova and causes the overall supernova 

^ I emission to be polarized. The motivation for the bipolar jet model is to resolve the 

O ! potential paradox set by an analysis which suggests strong spherical symmetry 

Q i for Type lb supernovae and the expectation that Type lb supernovae should 

r-| ! be highly polarized in the continuum up to perhaps ~ 4%. For the model we 

J^l have derived simple, non-relativistic, approximate analytic formulae for relative 

O I amount of scattered flux and polarization. The polarization formula qualitatively 

'^ I reproduces two main features observed for at least some supernovae: (1) the 

.. I rise to a polarization maximum and then decline with time; (2) the inverted 

.5^ I P Cygni polarization profiles of lines. The model is strictly axisymmetric and 

rS [ predicts a constant position angle of intrinsic polarization. Since position angle 

c^ I variation of intrinsic polarization is observed, the model cannot account for all 

the polarization of all supernovae. 

Subject headings: supernovae: general — polarization — radiative transfer 



1. Introduction 

1.1. Supernova Polarization in General 

Evidence has accumulated that net emission from core-collapse supernovae is often, per- 
haps usually (or maybe even nearly always), intrinsically polarized. So far the observations 
show a range in intrinsic continuum polarizations from ~ % up to perhaps 4 % with of 
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order 1% being typical (e.g., Wang et al. 1996, 2001; Wheeler & Benetti 2000; Leonard et 
al. 2000a,b; Leonard & Filippenko 2001; Leonard et al. 2001). (The ~ 4% polarization was 
actually obtained for a Type Ic supernova [Wang et al. 2001, but we assume it here as a 
tentative general upper limit on core-collapse supernova continuum polarization.) P Cygni 
line trough polarization features can exhibit higher polarization than the continuum: at 
the moment 4 % for Ca II line trough of SN 1987A seems to be the record (Jeffery 1991b, 
Fig. 16). The apparent polarization trends are that polarization increases with time from 
explosion for perhaps up to of order 200 days or more and that it increases with decreas- 
ing hydrogen/helium envelope mass of progenitor at the time of explosion. The sequence 
of decreasing hydrogen/helium envelope mass in core-collapse supernova types is II-P, II-L, 
Il-n, lib, lb, and Ic (Leonard et al. 2001). The order of II-L, Il-n, and lib is uncertain. 
For an explanation of the supernova classification scheme see reviews by, e.g., Filippenko 
(1997), Wheeler & Benetti (2000), and Cappellaro & Turatto (2001). (Type la supernovae 
[thermonuclear supernovae] seem to be usually much less polarized or unpolarized [Wang et 
al. 1996; Wheeler & Benetti 2000] although some detections have been reported [Wang et al. 
1997; Leonard et al. 2000b; Howell 2001]. We will not discuss Type la supernovae further in 
this paper.) 

The presence of intrinsic polarization shows some sort of asymmetry exists in the su- 
pernova ejecta or the near-supernova environment and that there is some sort polarizing 
scattering off electrons in the ejecta or, perhaps, dust well away from the actual ejecta. 
The continuum opacity of supernova atmospheres in the photospheric epoch in the optical 
and infrared (IR) is in fact dominated by electron scattering opacity (e.g.. Wagoner 1981; 
Harkness 1991), and thus the atmospheres themselves would yield a net polarization if they 
were asymmetric somehow. Ellipsoidal atmospheres models exploiting electron scattering 
have been extensively explored (e.g., Shapiro & Sutherland 1982; McCall 1984, 1985; Jeffery 
1991a; Hoflich 1991; and Howell 2001). Although, electron scattering has virtually no wave- 
length dependence for photons with energy much less than the electron rest mass energy 
(e.g., Davisson 1965, p. 51-52), strong line polarization features caused by the combination 
of electron scattering and line radiative transfer can occur in these models through a mech- 
anism investigated by Jeffery (1989). Another asymmetric supernova model invokes clumps 
in the ejecta as the cause of asymmetry (e.g., Chugai 1992). 

There is a model of supernova polarization that does not require an asymmetric super- 
nova. Instead it invokes a relatively remote scattering from a polarizing dust cloud (Wang 
& Wheeler 1996). Because the cloud is remote the scattered spectrum is added to the di- 
rect supernova spectrum of a later epoch: the scattered spectrum is time-delayed. Dust 
scattering has a relatively weak polarization wavelength dependence (e.g., Serkowski 1973), 
but the time-delay effect can give rise to line polarization features due to the variation of 



-3- 



the line profiles over the time delay. (In the bipolar jet model the line polarization features 
arise through a significant redshift of the jet-scattered flux from its original observer frame 
wavelength. We call this the redshift effect: see §§ 3.3-3.5 and 4). 

A considerable problem in analyzing intrinsic supernova polarization is contamination 
by interstellar polarization (ISP) caused by intervening dust along the line of sight from 
the supernova. That intrinsic polarization is present for almost all core-collapse supernovae 
is proven by the time dependence of the polarization (which would not happen for pure 
ISP) and the non-ISP wavelength dependence of polarization in particular across lines. (The 
ISP wavelength dependence is slowly varying and follows the empirical Serkowski law: e.g., 
Serkowski 1973.) Unfortunately, ISP can completely distort an intrinsic polarization spec- 
trum: e.g., polarization maxima can be turned into minima and vice versa. The correction 
for ISP must be based on the observed supernova data itself and, alas, on usually somewhat 
uncertain assumptions. Approaches to ISP correction are discussed by, e.g., Jeffery (1991b), 
Leonard et al. (2000a), and Wang et al. (2001). 



1.2. A Supernova Polarization Paradox and the Bipolar Jet Model 

Type lb supernovae (SNe lb) lack conspicuous hydrogen lines although most or all 
may exhibit weak hydrogen Balmer lines (e.g.. Branch et al. 2002, hereafter BBK). These 
supernovae may be essentially helium cores of a few solar masses that result from mass loss 
by massive stars (> 8 Mq on the main sequence) to a binary companion or in strong winds 
(e.g., Woosley & Weaver 1994, p. 122ff). The hydrogen remaining may be a few tenths of 
a solar mass or less which in the supernova ejecta is on the outside moving at velocities 
> ll,000km/s (BBK). (Type Ic supernovae are thought to be similar to SNe lb, but they 
have lost even more envelope and show no conspicuous optical helium lines although some 
helium remains in some cases at least [Filippenko et al. 1995; Clocchiatti et al. 1996]). 
The rise time to V maximum for low-envelope-mass core-collapse supernovae like SNe lb is 
of order 20 days as suggested by earliest observations (e.g., BBK) and theoretically (e.g., 
Shigeyama et al. 1994; Baron 2001). 

Recently, BBK have presented a parameterized LTE analysis of SN lb spectra. Out of 
a small, but perhaps representative, sample of 6 SNe lb BBK find that 5 show photospheric 
velocity curves (i.e., plots of the photospheric velocity evolution with time) that for the 
most part vary by < 20 % from a mean photospheric velocity curve for up to about the first 
80 days after explosion. Since photospheric velocity evolution is in fact a function of many 
properties of a supernova (e.g., mass, density distribution, composition, thermal state), the 
similar-evolution result strongly suggests that 5 out of 6 supernovae in this sample are quite 
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alike: tentatively the result suggests that most SNe lb could be quite alike. Moreover, the 
members of the subsample of 5 (hereafter the BBK subsample) are probably fairly spherically 
symmetric or else different viewing orientations would give different photospheric velocity 
evolutions even if the supernovae were identical. 

Recall now that SNe lb according to the cited trend of increasing polarization with de- 
creasing envelope mass and some observations (McCall 1985; Leonard et al. 2000b) should be 
highly polarized and have intrinsic continuum polarizations of order 1 % or perhaps some- 
times more. To obtain polarization of > 1 % from models with ellipsoidal density distri- 
butions requires for major-over-minor axis ratios > 1.2 % in the iso-density contours (e.g, 
Hoflich 1991; Hoflich et al. 2001). Such axis ratios suggest > 20 % variations in photo- 
spheric velocity curves for identical supernovae seen at the same epoch at different viewing 
orientations. It seems marginally unlikely that BBK subsample could exhibit the expected 
polarizations if ellipsoidal asymmetry is assumed. 

Unfortunately, only SN 1983N out of the BBK subsample has been measured for polar- 
ization. Spectropolarimetry for SN 1983N was briefly reported on by McCall (1985). (The 
SN 1983N spectropolarimetry data has never been published and now may not be recov- 
erable [McCall (2001)].) A polarization dip from 1.4% to 0.8% was found across an Fe II 
blended line flux emission near 4600 A near V maximum: the position angle of polarization 
stayed constant across the emission. The change in polarization across a line profile is a 
clear signature that intrinsic polarization was present, not just ISP (see § 1.1). The lack 
of variation in position angle suggests, in fact, that nearly all the polarization was intrinsic 
and thus that the intrinsic continuum polarization may have been as high as 1.4%. (Any 
significant change in the relative amounts of intrinsic polarization and ISP would tend to 
cause a shift in position angle, unless coincidentally the intrinsic and ISP position angles 
were aligned.) Thus, at least one of the BBK subsample is significantly polarized and it is 
the possible they all are. 

We face the potential paradox of fairly spherically symmetric supernovae yielding high 
polarization. To resolve this paradox, we propose a bipolar jet model for the polarization of 
core-collapse supernovae. In the model, bipolar jets are thrown out of the supernova at the 
time of explosion. If the jets are thrown out of the supernova core, it is plausible that they 
are enriched in radioactive ^^Ni. The explosion is such that these jets are well collimated 
and the bulk of the ejecta (which for short we will call the bulk supernova) is left fairly 
spherical: a contrived picture, but it is what we need for the resolution of the paradox. Since 
the bulk supernova is fairly spherically symmetric there is no viewing orientation variation 
in photospheric velocity provided one is not looking nearly down a jet axis. The jets stay 
sufficiently ionized (perhaps due to radioactive decay of ^^Ni and its daughter ^^Co) so that 
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their electron optical depth is high enough that light emitted by the photosphere of the bulk 
supernova and scattered by the jets into the line of sight is highly polarized. The photosphere 
and jet are sufficiently remote that the polarization position angle of the scattered light is 
nearly constant and perpendicular to the plane defined by the jet axis and the line of sight. 
(See the discussion of the polarization properties of electron scattering in § 3.3.) This means 
that the polarization of the scattered light is not completely canceled and may amount to 
a polarization of tens of percent. The polarized scattered light from the jets is only a small 
contribution to the net supernova fiux, but as we will show in § 4 that contribution can give 
the net fiux a significant overall polarization: e.g., > 0.5%. 

The existence of bipolar jets originating in core-collapse events has been much investi- 
gated recently and models have been calculated that show their effects (e.g, Hoflich et al. 
2001; MacFadyen et al. 2001; Zhang et al. 2003). The jets have been inserted into these 
models as parameterized entities. It is not clear if they arise naturally in core collapses. 
Thus, we are free to invoke bipolar jets with the properties we require without justifying 
them by the results of hydrodynamic calculations. 

In § 2 of this paper, we present the structure of our bipolar jet model. Section 3 describes 
the radiative transfer in the model and derives simple analytic formulae for relative amount of 
jet-scattered flux (i.e., the relative scattered flux) and the supernova polarization. We assume 
non-relativistic radiative transfer as a simplifying hypothesis. Because the jet is somewhat 
remote from the supernova, the scattered spectrum is added to the direct supernova spectrum 
of a later epoch: thus there is a time-delay effect that also occurred in the dust cloud model 
of Wang & Wheeler (1996). The time-delay effect could be included in our model, and may 
sometimes be important even for non-relativistic jets. We also allow for time evolution as 
the jets expand. Our simpliflcations allow us to capture the main effects of non-relativistic 
bipolar jets including time evolution. There are 5 independent basic parameters for the 
bipolar jet model: 6 (or alternatively /i), A^, Trad, An, and h. The Trad parameter can be 
reparameterized in terms of physical parameters for the jet (see § 2). In a calculation in 
which the time-delay effect is included, a reduced time, tred (see § 3.2), replaces Trad as one 
of the basic parameters. In a calculation of the time evolution of polarization, the parameter 
Trad.o^o replaces Trad (see § 2). These parameters are all explained in §§ 2, 3 and 4: see 
especially § 3.5. 

Section 4 reports some demonstrations calculations obtained using our analytic formulae. 
Conclusions and flnal discussion are given in § 5. 
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2. The Structure of the Model 

We imagine well-collimated, relatively narrow, conical, bipolar jets that emerged from 
the supernova core during the explosion event. The apex of the cones is at the supernova 
center, but we picture the jet cones as truncated at the inner and outer ends. The inner 
and outer ends we assume are sections of spherical shells concentric about the center of the 
ejecta. This assumption allows us to easily relate jet parameters and jet mass and kinetic 
energy as we do in this section and gives a simple mental picture for our formulation. In 
reflection phase of the jets (see § 3.3), we will, in fact, make the contradictory assumption 
that the inner ends are hemispheres with round sides toward the observer. This may be 
physically plausible and it permits an analytical treatment. The contradictory assumptions 
about the inner ends is not a concern. The structure of the real bipolar jets may be complex 
and we make assumptions that allow us to derive simple analytical formulae that we hope 
capture the main effect of such jets. 

The conical shape of the jets seems physically reasonable for material ejected from the 
supernova in an axisymmetric fashion as soon as the coasting phase (homologous phase: see 
below) has set in. We assume the jets are symmetric about the supernova center. The bulk 
supernova is assumed to be spherically symmetric. The bulk supernova and the jets con- 
stitute an axisymmetric system and thus the polarization position angle will be a constant. 
Given the nature of electron scattering (see § 3.3), this position angle will be perpendicular 
to the jet axis projected on the sky. (Position angle shifts of 90° are possible in general for 
axisymmetric systems, but will not occur for the bipolar jet model. See the discussion of the 
jet polarization in § 3.3.) 

The model as just sketched is not, of course, strictly applicable to jets in general. Still 
one can apply the model to jets in general with the understanding that the parameter values 
must be treated as characteristic values, not well-deflned values. Thus, the model can be 
applied to observations even if one is not assuming any precise structure for the jets which 
may for instance be more like clumps or fllaments. Simple variations on the model can also 
developed for analyzing actual observations. 

We consider only the homologous expansion epoch of the ejecta and assume that the 
jets are participating in the homologous expansion. In homologous expansion, the radial 
position Tj of any mass element i is given by 

n = Vit = Pict , (1) 

where Vi is the element's constant velocity, (3i is this velocity in units of c, and t is the time 
since explosion. (At explosion the supernova is usually effectively a point compared the size 
of the ejecta at observable epochs, and so initial radii are neglected.) Thus in the homologous 
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expansion epoch, velocity is a good comoving-frame coordinate for the ejecta. In this paper 
we will usually write velocity in units of c: i.e., we use /5's. 

The inner and outer radial velocities of the jets at their inner and outer ends are /5in and 
/9in + A/9rad, respectively: A/9rad is the radial velocity length of a jet. The velocities /3in and 
A/9rad are time-independent for the homologous expansion epoch, of course. The velocity of 
the supernova photosphere is /5ph: the photosphere is entirely in the bulk supernova. The 
photospheric velocity is time-dependent since the photosphere recedes into the ejecta as the 
density and opacity decline with expansion. We assume that Pi^ >> /3ph: i.e., the jets are 
well detached from the bulk supernova. This will allow us to approximate the photosphere 
as a point source. 

The angle of the jet axis from the line-of-sight axis in the direction of the observer is 6. 
We will often represent this angle by its cosine n = cos 6. We restrict 6 to the range [0, 90°], 
and thus fi to the range [0, 1]. The near jet to the observer is at fi and the far jet at — /i. 

The half-opening angle of a jet is A^. The fractional solid angle fifj. (i.e., solid angle 
divided by An) at the supernova center subtending a jet is given by 

f^fr(A^) = i[i - cos(A^)] ^ ^Ae' (i - 1a^^^ ^ ^Ae' , (2) 

where the second expression is the 5th order good approximation to the first and the third 
is the 3rd-order good approximation to the first expression. The 3rd order good expression 
(used with A^ in radians) is accurate to better than ~ 2.3 % for A^ < 30° and ~ 5.3 % for 

A^ < 45°. 

If we assume that the member of the BBK sample of 6 SNe lb that did not conform 
to the mean photospheric velocity curve was non-conforming because the line of sight from 
the observer to the supernova center passed through a jet, the estimate of the probability of 
the line of sight passing through a jet is 1/6. (We have no great confidence in this estimate 
because of the smallness of the sample and the fact that the bipolar jet model may not be 
right.) If we assume that all bipolar jet supernovae have the same A^, then the estimated 
probability implies that 2Qf^{A6) ^ 1/6 or A^ of order 35°. (The 2 factor is because our 
model has two jets.) We will not, however, take 35° as a fiducial half-opening angle A^: it 
is not at all a definitive result. The assumptions we will make about the scattering from the 
jet in § 3 are, in fact, better the smaller A^. 

We will assume we can treat the jets as geometrically narrow since this allows the picture 
of jet scattering we rely on in § 3.4. For this assumption to be really true 

A^An « A/?,ad or A^ « 6 , (3) 



where A^ is in radians, of course, and 
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The parameters b and /Sjn are chosen as two of the independent basic parameters of the 
model. The parameter A/^rad is a dependent parameter with this choice. We do study the 
behavior of the relative flux and polarization formulae (see § 3) for A^ values that are fairly 
large and cases with A^ > 6 in § 4 with the understanding that our assumed pictures of 
scattering in the jet (see § 3) are becoming poorer and must become invalid at some point 
as A^ grows large or b small. 

We will neglect occultation effects in our model. Occultation (partial or total) will occur, 
of course, for 6 sufficiently small: i.e., the near jet will occult the supernova photosphere and 
perhaps the far jet; the supernova photosphere will occult the far jet. Given our assumption 
that the jets are well detached (i.e., /^jn >> Pph), occultation should occur only for relatively 
small 6: in the limit of /9ph — > 0, occultation would only occur for 6 < A6. Our expression 
for polarization will in fact go to zero for 6 going to zero which is also what it should do if 
occultation were considered: the supernova and jets are axisymmetric about the line of sight 
when ^ = 0, and so cannot yield a net polarization whether or not occultation is considered. 
There may well be distinct effects on the line spectra caused by jet occultation. But those 
effects and the occultation situation in general are beyond the scope of this paper and we 
leave them for future consideration sine die. 

We assume the jets have a spatially constant electron density, density and ionization 
state. We also assume the ionization state is constant in time which means their electron 
density falls as t^^ due to homologous expansion. The jet electron density is 

' ma/ie ma/iefifr(A^)(47r/3) [{b + 1)3 - 1] f^lcH^ ' ^ ^ 

where p is mass density, ma = 1.66053886 x 10^^^ g is the atomic mass unit, M is the total 
mass of a jet, and /Xe is the mean atomic mass per electron which is determined from 

(e.g., Clayton 1983, p. 84): for ion i, Xj is mass fraction, Zi ionization stage, and Ai atomic 
mass. Assuming constant ionization state in time implies that He is constant in time too. 

The radial optical depth of a jet is 

Trad = neCTgA/^radCt (7) 
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where ae = 0.665245873 x 10~^^cm^ is the Thomson cross section for electron scattering and 
the final expression is in terms of fiducial values with r2fr(A) approximated by A^^/4 (which 
is accurate to better than ~ 2.3 % for A^ < 30° and ~ 5.3 % for A^ < 45° as mentioned 
above). We emphasize that r^ad decreases monotonically with time and becomes constant 
(in fact 0) only when time t goes to infinity. 

At present there is insufficient guidance for the fiducial values. The total mass ejected 
by typical SNe lb is perhaps only 5Mq or less (e.g., Hachisu et al. 1991), and one might guess 
the bulk spherical symmetry cannot be maintained unless a jet mass was of order tenths of a 
solar mass or less: hence we choose 0.1 Mq as the fiducial jet mass. The ejecta is expected to 
be helium enriched with intermediate mass elements, and perhaps some hydrogen: we have 
chosen the fiducial /ig to be that of singly ionized helium with the helium atomic mass set to 
4 exactly: thus the fiducial /ie = 4. The fiducial A^ = 20° ~ 1/3 radians and 6 = 1 are just a 
round numbers consistent with A^ << b (with A^ in radians of course) (i.e., consistent with 
geometrically narrow jets) that turn out to yield fairly high polarization (see § 4.2). The 
fiducial time t = 20 days was chosen since that is the order of the rise time to V maximum for 
low-envelope-mass core-collapse supernovae (see § 1.2). The fiducial An = 0.06 (i.e., inner 
end velocity ~ 18, 000 km/s) is about twice the photospheric velocity at V maximum and 
about three times the photospheric velocity at 30 days after V maximum according to the 
mean photospheric velocity curve of the BBK subsample. The choice of 0.06 as the fiducial 
An is consistent for V maximum epoch (marginally at least) and for later times with our 
assumption that An >> /3ph which allows us to approximate the photosphere as a point 
source. 

We note that with our fiducial values, the jets will have considerable kinetic energy. For 
one jet, a characteristic kinetic energy is 



1 . f M 

^char ~ -Mvl^, ^ 0.724 X 

char char ^O.IM© 



' An(l + fe/2) ^^ 

0.09 



foe , 



where Wchar = An(l + b/2)c is a characteristic velocity, foe stands for 10^^ ergs (the standard 
unit of supernova explosion kinetic energies), and we have used the fiducial values for M, 
An, and b (see eq. [7] above). Since the whole kinetic energy of an ordinary supernova 
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is of order 10^^ ergs (e.g., Woosley & Weaver (1994), p. 90), our fiducial values and t>char 
value require that the jets make large energy demands on the explosion mechanism. We 
are not concerned about this here. First, our fiducial values and fchar value were chosen 
without much guidance: perhaps less demanding values will turn out to be better in actual 
polarization spectrum analyses. Second, the actual energy in neutrinos available in core- 
collapse explosions is of order 10^^ ergs (e.g., Woosley & Weaver (1994), p. 90). It is thought 
that some of this available energy is transformed into the kinetic energy of the ejecta and 
is the main energy of the explosion. But how this is done exactly and how much energy is 
transformed are still very uncertain (e.g., Janka et al. 2004). Thus, we are plausibly free to 
invoke jet energies of order foes for our hypothetical jets. Hydrodynamic jet calculations also 
somewhat freely invoke jet energies (e.g, Hofiich et al. 2001; MacFadyen et al. 2001; Zhang 
et al. 2003). 

For our choice of fiducial values, we obtain jet radial optical depth 2.309 which is suf- 
ficient to produce fairly large polarization: see § 4.2. However, as equation (7) shows the 
optical depth decreases as 1/t^. Thus for increasing time with the other parameters fixed the 
optical depth falls and eventually polarization must decrease because few photons are being 
scattered and polarized. For example by day 100 (after the explosion) the radial optical is 
0.09 from equation (7). Such a low optical depth would usually give very low polarization 
compared to earlier times. But significant polarization for days 100 to 200 can be obtained 
from a different choice of parameters for the jets. The plausibility of these choices needs to 
be investigated. 

In an analysis where one does not want to specify the jet mass or mean atomic mass 
per electron /Xg, one can just specify 

Trad = Trad,0 ( y ) ' (9) 

where Trad.o is a fiducial radial optical depth at a fiducial time to. Actually there is only one 
parameter Trad.O'^O' t>ut splitting rrad,o^o ^^^'^ ^^'^ factors is an intelligible reparameterization: 
one might chose rrad,o = 1 in all cases and then to would become the single parameter. In 
this paper, we will often just use r^ad itself as a parameter and for time just a reduced time 
defined by 

tred = —;= with the inverse relation being Trad = -^— ■ (10) 

Note that when r^ad = 1, tj-cd = 1- The tj-cd = 1 epoch is important for jet polarization. Once 
the optical depth through a jet falls below ~ 1, the fraction of direct flux being scattered by 
the jet must start declining, and thus so must the jet polarization (see § 4.2). 
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Using equations (7) and (10) we find the relation of time to reduced time to be 

h 1 



t 



iredy Trad, 0^0 — ^ 



3 a 



M 



rcdi 



47rm,/iefifr(A^)[(6+l)3-l]A„c 
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0.1 Mo 
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(12) 



where the final expression is in terms of our fiducial values for the parameters and we have 
approximated i7fr(A) by A^^/4 again. If the fiducial values are used, we see that tred = 1 
and Trad = 1 occur at day 30.39 after the explosion. 

Another quantity we will need is the radial velocity, which we will call /3scat, which 
locates the place in the jet that is the mean optical depth into the jet that photons radially 
incident on the jet reach for their first scattering under the assumption that they do scatter 
in the jet. This mean optical depth Tmcan is given by 

Trad 



' TT 



1- 



exp(ri.ad) - 1 



-T"rad I 1 



^rad 

~6~ 



r; 



rad 

360 



in general; 
for Trad >> 1; 

for Trad < 1- 



(13) 



At Trad = 0, Tmean = 0; Tmean then grows mouotouically with Trad initially with slope 1/2; but 
the slope decreases monotonically with r^ad and rmcan goes asymptotically to 1 as r^ad -^ oo: 

T-mean(T-rad = 1) ~ 0.418, T^cmi{rvaA = 3) ^ 0.843, and rmean(T-rad = 5) ^ 0.966. The Ouly 

stationary point of Tmean is a maximum at r^ad = oo. In supernova time evolution, Tmean is 
nearly a constant 1 at early times when r^ad is large, and then decreases as r^ad decreases; 
Tmean is decreasing roughly linearly with r^-ad when r^ad falls to of order 1; then Tmean vanishes 
linearly with r^ad as time goes to infinity. 



Since the electron density is assumed uniform, /5scat is given by 



Pscat — Pin + Aprad" 



/ rr 



AA 



rad 



,-1 



&A. 



,-1 



/ rr 



'rad 



(14) 



'rad \ 'rad 

We take /3scat to be the characteristic comoving-frame radius for scattering in the jet. At 
late times when r^ad ^ 0, Aeat -^ An + AArad/2. 

We will also need the optical depth perpendicular to the jet axis from the jet axis to the 
jet surface at Ascat- We call this the off-radial optical depth Tos- The approximate expression 
for Toflf (with A6' in radians) is 

'rad 



Toff = A^A 



scat 



AA 



rad 



M{^-f 



+ T-n 



(15) 



-12- 



For high polarization, it is good to have r^ad large for lots of polarizing scattering and Tos 
small for low depolarizing multiple scattering (see § 3.4). This means A9/b << 1 is a 
desirable condition for high polarization: we have, in fact, already indicated our need for 
this condition (which implies a geometrically narrow jet) above (see eq. [3]). In another 
perspective. Trad large and Tqs small together imply a geometrically narrow jet since we are 
assuming constant electron density in the jet. 

Mathematically, r^s always decreases with decreasing r^ad or increasing time. Its behav- 
ior is summarized in the following expressions: 



(t^ + T"mean j CXactly; 



^rad 
''~rad 



Toff ^ A^ I 



for early times when r^d >> 1 and r^ad » b; 
for all times if 6 << 2; 



r 



Trad ( T + 7T 1 fo^ ^^^^ timCS wheU Trad S, 1 



L for intermediate times if 6 >> 2; 

Trad/^ << 1, and Trad >> 1- 

Note the following also. If 6 < 2, then Trad/^ > Tmcan always, except at time going to infinity 
where both Trad and Tmcan are zero. If 6 = 2, then r^ad/b and Tmcan asymptotically approach 
each other as time goes to infinity. If 6 > 2, then there is a phase where r^eud/b < Tmcan, that 
ends only when time goes to infinity. 

The time evolution of r^ad and Toff control the time evolution of bipolar jet model polar- 
ization as we will see from the analytic polarization formula equation (48) and equation (42) 
in § 3.5. An interesting prediction of the bipolar jet model emerges from this control is 
there are two kinds of phases in which the continuum polarization can plateau for some time 
interval if we assume that no significant time dependence affects the continuum polarization 
through the redshift effect (see § 3.3). The first kind occurs at early times going back to 
time zero when Trad and Tofr are both effectively infinite for a finite time: effectively infi- 
nite since they enter the function controlling polarization evolution through the functions 
1 — exp(— Trad) and exp(— Tofj) (see § 3.5, eq. [42] and [48]) and make those functions nearly 
exactly 1 and 0, respectively for a finite time. This kind of polarization plateau phase is 
demonstrated in the calculations in § 4.2: see the early-time, non-infinite b polarization 
evolution curves in Figures 1, 2, and 3 in § 4.2. 

The second kind of polarization plateau phase occurs at intermediate times for b » 2 
or back to time zero for b = oo. In this case, r^ad is effectively infinite in the functions 
1 — exp(— Tj-ad) and recall making both these functions nearly 1 for a finite time, but, on the 
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other hand, Trad/^ << Tmcan ~ 1 for a finite time (which situation requires 6 >> 2) and this 
makes Tog ~ A^ (with A^ in radians, of course) for a finite time. In other words, one requires 
Trad/b « I while Trad >> 1 for the second kind of polarization plateau. The second of kind 
of polarization plateau phase is also demonstrated in § 4.2: see the polarization evolution 
curves with 6 > 100 in Figures 1 and 3. 

That the polarization plateau phases can exist is a mathematical prediction of the 
bipolar jet model. The first kind of plateau phase is unavoidable if our assumptions are 
adequate and should be observable if a supernova with jets is observed early enough. The 
second kind can only exist if 6 >> 2: we show calculations without the second kind of 
polarization plateau in in Figures 2 and 3 in § 4.2. The second kind of polarization plateau 
is probably physically unrealizable for cases of significant polarization as we discuss in § 4.2. 



3. The Radiative Transfer of the Model 

In the derivation of the bipolar jet model formulae for scattered flux, relative scattered 
flux, and polarization, we make many simplifying assumptions. These assumptions are 
needed to capture the main effect of the model on flux and polarization analytically and 
tractably. One main simplifying assumption is that the only signiflcant opacity in a jet is 
electron scattering (or Thomson) opacity. Thus we neglect the scattering, absorbing, and 
depolarizing effects of other opacities of which line opacity is probably the most important. 
Another consequence of assuming only electron scattering opacity is that no photons are 
destroyed in the jet. We also assume none are created. 

In § 3.1, we will make the point source assumption for the spherically symmetric bulk 
supernova. Subsection 3.2, gives a simplifled treatment of the time-delay effect which could 
be incorporated into the bipolar jet model. In the following three subsections, we derive 
analytic formulae for scattered flux, scattered relative flux, and polarization spectra that 
account for the three phases of jet polarization. The flrst phase, which we call the reflection 
phase, is when the jet is optically thick in all directions (see § 3.3). The second phase, which 
we call the radial phase, is when jet is optically thick in the radial direction, but optically 
thin in the direction perpendicular to the jet axis (see 3.4). The third phase, which we call 
the optically-thin phase, is when the jet is optically thin in all directions (see § 3.4). In § 3.5, 
we give general expressions that are for all phases. 
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3.1. The Point Source Approximation for the Supernova 

Given that the direct luminosity from the supernova (i.e., the luminosity without the 
jet contribution) at rest-frame wavelength A is Ldir(A), the direct flux at Earth from the 
supernova (i.e., the flux without the jet contribution) is 

b-W = ^ . (17) 

where D is the Earth-supernova-center distance. 

The fraction of Ldii.(A) that enters the inner end of a jet is just Qf^{A6) (see eq. [2] in § 2). 
This is true regardless of the size of the supernova photosphere (provided /9ph < An): it is a 
consequence of the assumed spherical symmetry of the spherical bulk of the supernova. We 
will assume that Ld[T.{X)Q{j.{A6) is all the luminosity that passes through a jet and that all the 
beams from the bulk supernova are radial: i.e., all the beams from the bulk supernova that 
enter a jet are radial and they enter at the inner end and they leave at the outer end. Thus we 
are making the point-source approximation for the photosphere. This approximation should 
be reasonably good as long as Pi^ » /3ph as we have already assumed (see § 2). There are, 
of course, two errors associated with the approximation caused by non-radial beams from 
a flnite photosphere. First, there will be non-radial beams entering the inner end of a jet 
(these are counted in Ld[^{X)Qfj.{A6) , of course) and these beams will not always traverse 
the whole radial extend of the jet. Second, there will be non-radial beams that enter the 
sides of a jet (these are not counted in Ldir(A)r2fi.(A^), of course) and these beams will only 
traverse some part of the radial extent of the jet. The flrst error leads to an overestimate of 
and the second error, to an underestimate of the direct-flux-jet interaction when we make 
the point-source approximation. The two errors are somewhat compensating: this somewhat 
strengthens point-source approximation. The point-source approximation is also somewhat 
strengthened by the fact that the photospheric emission will be somewhat limb-darkened: 
thus the photosphere will be more point-like than assuming that angle-independent speciflc 
intensity beams emerge from it. Thus, our photosphere is more point-like than a non-zero 
/9ph suggests at flrst. 



3.2. The Time-Delay EflPect 

The time-delay effect is caused by the supernova evolving during the time that flux 
travels from the photosphere to the jets and during the radiative transfer through the jets. 
In this section we will develop some formalism to treat the time-delay effect in the bipolar 
jet model, but we will not actually incorporate this formalism into our formulae in §§ 3.3, 
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3.4, and 3.5 and our calculations in § 4. The variations on those formulae that allow for the 
time-delay effect can easily be derived if needed, but are a bit cumbersome in appearance. 

In general, time-delay effect will be complex. So we will make some simplifying as- 
sumptions. We will assume again that the bulk supernova can be treated as a point at the 
center of the system and we will assume that the jet can be treated as a point at velocity 
Acat (see § 2, eq. [14]) for time-delay effect purposes: this eliminates the need to consider 
the time delays while flux transfers through the jets. Now /Sgcat itself is time dependent. As 
a simplification we assume it can be evaluated at a fiducial time t: t and all other times 
below are counted from the time of the supernova explosion: i.e., since time zero. This last 
assumption is actually rather poor since the time-delay effect will cause us to evaluate jet 
optical depth parameters at times different than t (as we will see below) and those optical 
depth parameters determine /Sgcat- But to get a sense of the time-delay effect, it is adequate. 

For each time of observation of flux and polarization spectra, we will in general have 
to consider 5 supernova evolution times. But we will compact the expressions by using an 
upper case label for the near jet and a lower case label for the far jet. Let tjet± be the time 
that flux arrives at a jet from the bulk supernova. Let t'-^_^ be the time that same flux 
started out from the bulk supernova. The direct flux from the bulk supernova that starts at 
flducial time t arrives at the Earth at the same time as the flux from the jet that starts at 
time tjct±: as mentioned above, we assume no time delay for the flux inside the jet. From 
geometry and homologous expansion, it follows that 

^jet± = t^eti H — and tjct± = t ± — , (18) 

c c 

where /x (as specifled in § 2) is the cosine of the angle to the near jet: we conflne fi to the 
range [0, 1]. A little algebra leads to 



^ Pscat 



t']et± = t-, |r^ and tjet± = t- — . (19) 



From the expressions of equations (18) and (19), we have at once the inequalities 

tj'et+ < t < ^jct+ and tj',t- < ijct- < t . (20) 

The inequalities show the time ordering of the events. Note that tjct± = t only when /U/3scat = 
and t'^ct± = ^jetit only when /5scat = 0. Note also that the ratios of the time parameters 
are not constant since /Sgcat is itself time dependent. We have recall simplifled this time 
dependence by assuming that we can evaluate /Sgcat at the flducial time t itself. 

Equation (19) is the probably the most useful time-delay effect result since flducial time 
t is the natural time for any calculation: it is the time that characterizes the bulk of the 
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flux that ones sees: sees D/c later than t, of course: recall D is the Earth-supernova-center 
distance. 

As we will see in §§ 3.3, 3.4, and 3.5, polarization is high for /x ^ 0. For example, for a 
high polarizing jet one might have /i = and /Sgcat = 0.1: this implies t'-^ct± = ^■'dt. Thus, with 
these parameters the observer detects a mixture of direct flux and scattered flux where the 
scattered flux originated in the bulk supernova at at time 10 % earlier than the time of origin 
of the direct flux. Since considerable evolution in flux level and spectra can occur in 10 % of 
the time since explosion, the time-delay effect could be signiflcant. Thus, in general it would 
be useful to have a complete enough time coverage of observed flux spectra so that one can 
use the observed spectra from earlier times to construct synthetic polarization spectra using 
the bipolar jet model for later times. 

Equation (19) can be used with reduced times (see § 2) replacing ordinary times: i.e., t 
becomes tred) ^jct± becomes trcdjeti, and tjg^_|_ becomes t^^^jg^.^. Using tred as the flducial time 
and equation (19), and equations (10), (13), and (15) of § 2, one can calculate the needed 
time-delay affected optical depths for trcd- Also using trcd and equation (19), one can calculate 
the times at which the flux from the jets started out from the supernova (i.e., t'^^^ jet±)' ^"^^ 
thus when calculating the net flux one can add the direct and jet fluxes corrected for the 
time delay effect. In fact, t^cd (which replaces t) can be used as of the 5 independent basic 
parameters of the bipolar jet model when the time-delay effect is included. If one neglects 
the time-delay effect. Trad is the same for both jets in a calculation and r^ad can be used 
instead of t^d as one of the 5 independent basic parameters. 

We will not explore the consequences of the time-delay effect in detail in this paper. 
But we will briefly discuss them here. As a heuristic formula (suggested by eq. [48] in § 3.5), 
let polarization as a function of t be given for the nonce by 

p(t)~Po(t)n(t), (21) 

where Po{t) is polarization in the absence of time delay and 

n(t) = ^-^ , (22) 

where /dir is the direct flux from the bulk supernova and tj^^ is some appropriate average 
of tjg^^ and tjgt„. As shown by the demonstration results in § 4, there will be a time of 
maximum (continuum) polarization. These demonstration calculations were done neglecting 
the time-delay effect, and so their polarization corresponds to what we have labeled for the 
nonce Po{t). Clearly, if n(t) is a rising/falling function at the time of the Po{t) maximum 
(which is a stationary point recall), then the P{t) maximum will be shifted to a later /earlier 
time than that of Po(^) maximum. 
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Now it is difficult to predict the temporal behavior of n(t) in general since it depends 
on the shape of the supernova light curve and on the how tL varies with time. But near 
maximum light, n(t) should be a rising function since the fdirit'jct) ^^^^ t)e rising and /dir(^) 
should be nearly stationary. So if the Po{t) maximum is at about maximum light, the 
time-delay effect will tend to give the P{t) maximum at a later time. 

At late times (i.e., probably of order 100 days or more after explosion), /3scat will approach 
a constant A/9rad(^"^ + 1/2) as Tj-ad -^ (see eq. [13] and [14] in § 2), and thus the ratio 
tjgj/t will approach a constant too (see eq. [19]): a constant that is less than 1, of course 
(see eq. [20]). At such times the light curve of the supernova is falling exponentially or 
quasi-exponentially (e.g., Jeffery 1999), and so n(t) ?s exp[(l — i'^^^/t)t/te] which is growing 
exponentially: te is the assumed e-folding time of the light curve and the ratio tLt/^ i^ a 
constant less than or equal to 1 (see eq. [19] and [20]). Thus if the Po{t) maximum occurs 
in this late phase, the P{t) maximum will tend to occur at a later time again. Actually, 
because ejecta is tending to be optically thin in the exponential/quasi-exponential phase, 
one does not expect the Po{t) maximum to occur then. 

In the analytic formulae we derive in the remainder of this section and in the calculations 
reported in § 4, we will, as mentioned above, neglect the time-delay effect. This means we 
approximate all the time parameters by one value, t or trcd- When only one time parameter is 
needed then, as mentioned above, r^-ad can be used instead of trcd as one of the 5 independent 
basic parameters. 



3.3. Reflection Phase 

In the reflection phase, the jet is optically thick and it acts as a diffuse reflector of 
incident light. Recall that for simplicity we have assumed that the jets have uniform electron 
density. In fact they are likely to increase in electron density as one moves inward. Thus 
photons striking the optically thick jet that penetrate somewhat deeply in optical depth will 
scatter a few times and then tend to diffuse outward again and be emitted from the jet. The 
behavior of these deep-penetration photons is effectively uninfluenced by their history prior 
to entering the jet. The deep-penetration photons will have an emergent angular distribution 
that is rather isotropic over solid angle 2it (i.e., the half of the sky not filled by the jet) and 
is mostly unpolarized — individual photons will be polarized, but their net flux mostly not. 
Some anisotropy and net flux polarization will occur due to diffusive photon scattering in an 
inhomogeneous medium (e.g., Chandrasekhar [1960, p. 248] for plane-parallel atmospheres; 
Cassinelli & Hummer [1971] for spherically symmetric atmospheres), but this is likely to be a 
relatively minor effect and too complicated for us to contemplate in the simple picture we are 
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developing. There is also a class of photons that scatter more than once, but not sufficiently 
to lose all trace of their history before they entered the jet. Their angular distribution and 
polarization state are still somewhat determined by their initial direction of incidence on jet: 
but this non-random redistribution and polarization is likely to have a relatively small effect. 
We will assume that photons that scatter more than once in the jet emerge unpolarized and 
are isotropically distributed over 27r. 

The photons that scatter once and escape the jet will have an angular distribution 
and polarization determined by the Rayleigh phase matrix (e.g., see Chandrasekhar 1960, 
p. 37). We will estimate the fraction q of photons in this category in the following way. 
The probability that a photon penetrates a jet along some typical beam path to optical 
depth r without scattering and then scatters in optical depth dr is exp(— r)(ir. Since we 
are considering a very optically thick phase, we can regard the local jet surface as a plane- 
parallel atmosphere. Thus ~ 1/2 the scattered photons are aimed at a path toward the 
"plane-parallel" surface and ~ 1/2 are directed inward and must scatter more than one 
time before escaping: we assume these fractions are both exactly 1/2. This half-inward- 
half-outward division among scattered photons is true even with the anisotropic Rayleigh 
distribution since this distribution is axisymmetric about the incident direction and forward- 
back symmetric: any plane through a polar plot of the distribution slices it into two regions 
with the same integrated emission. We will assume that the average optical depth distance 
to escape the jet for a photon that has just scattered for the ffist time is equal to r: note 
the penetrating and escaping directions are not the same; we are just assuming the optical 
depth along them is the same in some kind of average sense. Thus the probability of this 
photon escaping without scattering again is estimated to be exp(— r). If we multiply the 
relevant factors together, we find the estimate of the probability of a photon incident on the 
jet penetrating to optical depth region dr and then escaping the jet is (1/2) exp(— 2r) rfr. 
The estimated q fraction is the integral of (1/2) exp(— 2r) dr over all r: 

q^-j exp(-2r)rfr = -, (23) 

where the integral can be extended to infinite optical depth because we are assuming the jet 
is optically thick. If one had a well specified jet structure, one could calculate more exactly 
the average fraction of once-scattered photons. That calculation is too elaborate for our 
crude picture of the jet. In this paper, we adopt 1/4 as the fixed value for q. 

To treat polarization we must label polarization axes for the scattered light. We will 
label these axes (which are both perpendicular to the scattered light's propagation direction) 
X and y: the x axis is perpendicular to the plane of scattering and the y axis is aligned with 
it. The plane of scattering (or scattering plane) is the plane defined by incident and scattered 
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beams. The Stokes Q and U parameters are needed to evaluate the polarization: see see 
Chandrasekhar (1960, p. 25) for a detailed presentation of the Stokes parameters. The 
Stokes Q parameter for the scattered flux in our setup is |/-flux minus the x-flux following 
conventions of Chandrasekhar (see Chandrasekhar 1960, p. 27). 

Because of the finite size of the jet there will be a range of scattering planes for light 
scattered in the observer direction. In our mental picture, we will approximate this range of 
planes by the plane defined by the jet axis and the scattering direction to the observer: we 
will call this plane the fiducial plane. Just using the fiducial plane alone would overestimate 
the polarizing effect since the polarization-canceling effect from a range of planes is neglected. 
We will make a crude correction factor for the Stokes Q formula to account for the range of 
scattering planes. 

Consider a circular arc of material in the jet with arc center at the center of the supernova 
in the plane defined by the jet axis and the line perpendicular to the line of sight. The arc is 
meant to approximate the layer of scattering material in the jet at /Sscat- The angle at the jet 
center subtended by the arc is naturally 2/\9. Measure the angle for any point along the arc 
by C- C = si't the jet axis and is positive counterclockwise. Suppose again that the supernova 
photosphere is a point. Emission from the supernova center is scattered off electrons in the 
arc toward the observer. Assume the incident specific intensity on the arc, the amount of 
scattering, and the polarizing effect are constant along the arc. The local Stokes axes at each 
point along the arc are labeled x' (tangent to the arc) and y' (perpendicular to the arc). The 
unadorned x and y axis are fiducial axes for the fiducial plane. The local Stokes parameters 
are dQ' and dU': dQ' < from the Rayleigh phase matrix (e.g., see Chandrasekhar 1960, 
p. 37) and dU' = by symmetry: we are assuming the incident light is unpolarized. The local 
position angle of polarization is aligned with the x' axis: see the Rayleigh phase function 
(eq. [28]) and accompanying discussion below. If we could legitimately just integrate up the 
local dQ' without transforming to the fiducial x-y axes, then the net Q' for the arc emission 
would be proportional to 2A6. More correctly we should transform dQ' and dU' to the 
fiducial x-y axes and then integrate. The transformations for a rotation ( clockwise are 

dQ = dQ' cos(2C) and dU = -dQ' sin(2C) (24) 

(Chandrasekhar 1960, p. 34). 

The integration for net U gives zero because the integration interval is even: from 
a more obvious perspective, the axisymmetry of the jet requires the net U for the fiducial 
plane Stokes axes to be zero. Because U = 0, the position angle of polarization is determined 
only by the net Q. Given the nature of electron scattering (which we discuss below), the 
polarization is aligned perpendicular to the fiducial plane. 
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The net Q corrected for the variation of the polarization axes along the arc in terms of 
the net Q' neglecting this variation is given by 

Thus to crudely correct the Stokes Q parameter in our formulae developed assuming only 
fiducial plane scattering, we will multiply it by the correction factor 

e(A^) = ^^^^ - 1 - 1(2A^)^ + -^^{2Aer , (26) 

where the second expression is the 5th order good approximation in small 2A6 (in radians) 
to the first. Since the 1st order term in the correction factor vanishes, we would have 1st 
order good expression for Q without using the correction factor. The correction factor is, of 
course, only a crude correction since the jet is not a single arc at a single inclination to the 
line of sight. But the ^(A^) factor does crudely have the right effect and the exactly right 
limiting behavior. In the vanishingly-narrow-jet limit A^ -^ and ^(A^) ^ 1 as it should 
since there is no range of scattering planes in this case. In the limit where the jets cover 
the sky spherically symmetrically A^ -^ it/ 2 and .^(A^) -^ 0: this is also correct since the 
symmetry would cancel all polarization. Between the limits, C,{A.6) decreases monotonically 
with A^ with the only stationary point being a maximum at A^ = 0. At worst, ^(A^) is an 
adequate interpolation formula for a correction factor; it should be somewhat better than 
that though. 

A similar Stokes Q parameter correction factor could perhaps be devised to account for 
the finite size of the supernova photosphere. But it would complicate our mental picture 
and, given the crudity of the model, seems of little value. 

The angular redistribution function for electron scattering of unpolarized light is the 
normalized Rayleigh phase function: 

fi'norm(,A''gcnJ ~j 9\f^gcn) ■ [■^ ' ) 

g{figcn) is the (conventional) Rayleigh phase function given by 

3 

^(/igcn) = 4(1 + /^gcn) , (28) 

where Hgen is the cosine of a general scattering angle (i.e., the angle between incident and 
scattered beam directions) (e.g., Chandrasekhar 1960, p. 35). The net direct supernova light 
incident on the jets is unpolarized. Individual beams are polarized since the supernova at- 
mosphere is electron-scattering dominated, but spherical symmetry of the bulk supernova 
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makes the net polarization zero. In our notation the 1 in the 1 + yUg^^ factor determines 
the intensity perpendicular to the scattering plane (i.e., the x-flux) and the /Xgen, the inten- 
sity parallel to the scattering plane (i.e., the y-Hux): this follows from the Rayleigh phase 
matrix (e.g., Chandrasekhar 1960, p. 37). Note that forward and backward scattering are 
unpolarized since /Xgg^ = 1 in those cases and the two scattered flux components are equal. 
Scattering through 90° is gives the maximum polarization of 100 %. 

The Rayleigh phase function is adequate to describe photon redistribution for us for 
two reasons. First, we can assume electron scattering is coherent (i.e., does not change 
wavelength) in the jet frame. The coherent scattering approximation is adequate since 
the wavelength part of the photon redistribution (which affects the angular redistribution in 
general) on electron scattering for electrons obeying a Maxwellian distribution (Mihalas 1978, 
p. 420) averages away for continuum radiative transfer which is all we are considering in the 
jet. Second, the flux in the jet can be approximated as unpolarized in an average sense before 
a last scattering which may or may not be polarizing. For the reflection phase considered 
in this section the last scattering in the jet is only polarizing on average if it is also the first 
scattering as we argued above. For the radial and optically-thin phases of the jet (see § 3.4), 
the photons that are only scattering radially inward and outward (i.e., radially directed 
photons) are approximately unpolarized since there scattering is approximately backward 
and forward scattering. The last scattering is polarizing if it is of radially-directed, and so 
unpolarized, photon. If the last scattering of non-radially-directed photon is not polarizing 
on average as we will argue in § 3.4. 

We assume here and in all our developments that all scattering toward the Earth that 
is polarized (on average) can be treated as being from the jet axis to the direction toward 
the observer, and so the only cosines of scattering angles for this flux are fi and — /i: i.e., 
the cosines for the near and far jets from the line of sight in the direction to the observer 
(see § 2). (For the reflection phase this approximation is justified by our assumption of A^ 
small. We justify it for the radial and optically-thin phases in § 3.4.) The approximation is, 
of course, somewhat gross. The unpolarized (on average) scattered fiux toward the observer 
we approximate as coming from an isotropic angular distribution, and so the cosine of the 
scattering angle is unneeded and unspecified. 

Although we treat the electron scattering as coherent in the jet frame, the supernova 
fiux that scatters off the jet and is emitted toward the observer at rest-frame wavelength A 
is redshifted from the rest-frame wavelength A' that it had on emission from the supernova 
photosphere. We call this shifting the redshift effect. It is caused by macroscopic-velocity- 
field Doppler shifts. The approximate relation between the two wavelengths from the non- 
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relativistic Doppler formula is 

A = A'±[l+/5scat(lT/i)] , (29) 

where the upper and lower cases are for scattering off the near and far jets, respectively, from 
the observer: recall fi is the cosine of the angle of axis of the jets measured from the line of 
sight axis in the direction of the observer. The /Sgcat velocity was defined in § 2 (see eq. [14]) 
as the characteristic radial velocity of jet scattering: for the reflection phase /Sgcat ~ An since 
the jet is very optically thick, and so Tmcan/Trad << 1- There actually will be somewhat 
more redshift than given by equation (29) caused by multiple scattering in the homologously 
expanding medium of jet. In the reflection phase those scatterings are close together in 
velocity space, and so probably cause negligible extra Doppler shift; in later phases such 
shifts may be significant, but we will neglect them for simplicity. We actually will use the 
approximate inverse of equation (29) in our formalism below: 

A'± = A[l-/3seat(lT/i)] , (30) 

which is an inverse formula good to 1st order in Acat- 

To calculate the net jet-scattered flux and the net Stokes Q parameter for this flux in the 
reflection phase we will, as mentioned at the beginning of § 2, assume the inner ends of the 
jets as hemispheres with round sides toward the center of the supernova. This assumption 
may be physically plausible and it allows an analytical treatment. But it does contradict 
our basic picture of the jet ends as discussed at the beginning of § 2. 

Let the cross sectional area of the hemispheres be A. Toward the observer the reflecting 
face of the near jet hemisphere will be crescent-shaped and that of the far jet hemisphere 
will be gibbous. The faces of the hemispheres on sky have projected areas 

^(lT/i)A, (31) 

where the upper and lower cases are for the near and far jets, respectively. Note that the 
summed projected area of the jets is just A. 

Approximating the supernova photosphere as a distant point source, the hemispheres are 

wholly illuminated by the direct supernova flux. The direct supernova luminosity (emitted 

at wavelength A'_|_) absorbed per unit projected area (i.e., projected area as seen from the 

supernova center) by a jet is 

Ldi.(A-^)^fr(Ag) 

A • ^^^^ 

An approximate expression for the reflected specific intensity emitted by a hemisphere in 
the direction toward the observer is 

^^jr4 [^(/^)?+(l-?)] • (33) 
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where the g{fi)q term accounts for the once-scattered, polarized photons, and the (1 — q) for 
the muhiply scattered, unpolarized photons. Recall we assumed that we can approximate 
all cosines of scattering angles for the once-scattered photons by the fj, and — /x. We only 
need the /x in equation (33) for both jets since the phase factor formula only depends on the 
square of the scattering angle cosine (see eq. [28]). The denominator 2nA in equation (33) 
follows from the following considerations. If the energy absorbed by one jet from the direct 
flux from the bulk supernova were emitted by specific intensity beams that were constant 
with angle of emission and surface location from a sphere of radius r, then the denominator 
would be 47r^r^, where the second vr factor accounts for the difference between flux (using 
in the term in its formal radiative transfer sense) and specific intensity. In our case, all the 
absorbed direct flux energy gets radiated out of a hemisphere and not a sphere, and so the 
numerator is set to 27r^r^ or 2ttA taking r to be the radius of the hemisphere. 

Equation (33) is approximate because we have assumed a remote point photosphere and 
because of the approximate nature of the q factor. It is also approximate because we assume 
that the energy collected from the bulk supernova is spread evenly over the hemispherical 
surface. Actually the edges of the surface get less energy per surface area than the center 
because of the hemispherical curvature: the reflection specific intensity should decrease as 
one moves from center to edges. However, the overall reflected flux works out nearly correctly 
as we show below probably because of the compensation we get when summing up all the 
flux from both jet hemispheres emitted toward the observer. 

Now the flux at Earth from an object emitting specific intensity I is f = j I dVt = 
J I dA/D"^, where dfl is the differential bit of solid angle at Earth that the object subtends 
and dA is a differential bit of the projected area of the object on the sky (e.g., Mihalas, 
p. 11). Given this fact and making use of equations (17), (30), (31), and (33), it follows that 
the total jet flux measured at Earth at the rest-frame wavelength A for the reflection phase 
is 

/reflection (A) = [/dir (A'^) (1 " /i) + /dir(A'.) (1 + /i)] fifr (A^) 

[g{^)q+{l-q)] . (34) 

We note that if A'_,_ and A'_ are both set to A and q to zero, then 

_ WA)fMA^ 

/reflection I A j — ^ T — Jy^ • \O0} 

This last result is the flux at the Earth if all the luminosity captured by the jets were 
radiated isotropically. That we get this result verifies that we are accounting for the energy 
approximately correctly in equation (33) for the reflected specific intensity. 



-24- 



Froin equation (34) and making use of our understanding of equation (28) for the 
Rayleigh phase function, we find that the net Stokes Q parameter for the refiected fiux 
is 

Qrcflcction(A) = [/dir(A'J (1 - /i) + /dir(A'_) (1 + /i)] f]fr(A^) 

(^) (/i' - 1) q^^O) , (36) 

where we have introduced the correction factor ^(A6') (see eq. [26]) that accounts for effect 
of the range of scattering planes on the Stokes Q parameter. Note the {1 — q) term in 
equation (34) accounts for unpolarized flux, and thus the x- and y-fluxes it governs are equal 
and it cancels out of ^reflection- 



3.4. Radial and Optically Thin Phases 

In the ideal picture of the radial phase, a jet is optically thick in the radial direction, and 
so no photons can escape out the outer end of the jet. On the other hand, we assume the jet is 
optically thin perpendicular to its axis. With this setup, radially-outward scattered photons 
remain trapped in the jet as do the radially- inward scattered ones (except at inner end of the 
jet) and off-radially scattered photons are assumed to escape without scattering again. Recall 
that forward and backward scattering are unpolarizing as determined by the Rayleigh phase 
matrix (see § 3.3), and thus the radially scattered photons remain unpolarized. As mentioned 
in § 3.3, we assume that all scattering toward the Earth can be treated as happening along 
the jet axis, and so the cosines of scattering angles are again /i and — /i. Thus there is no 
averaging over a range of scattering cosines. We again will include our correction factor 
^(A6') (see § 3.3, eq. [26]) to account for the range of scattering planes. Recall the scattering 
plane deflned by the jet axis and the line of sight to the observer is the flducial plane. 

Given aforesaid picture, all the photons that escape the jet are distributed in angle 
and polarized as if they had undergone a single electron scattering event (i.e., their last 
scattering which was off-radial or radially-inward) — except that in our picture there is no 
flux scattered radially-outward and perhaps less scattered radially- inward since some photons 
should random walk far out in the radial direction. The exception tends to cause more 
polarized emission than predicted by the Rayleigh phase matrix for a single scattering. 
Trying to invent a correction factor for this effect that is of the right size seems futile given 
the crudity of the model. 

To handle the phase where the jet becomes optically thin radially (i.e., the optically- 
thin phase), we include the factor 1 — exp(— Trad) in our expressions for flux and Stokes Q 
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parameter: 1 — exp(— Trad) is the fraction of photons that scatter in the jet at least once: it 
goes to zero as Trad goes to zero. Thus the expressions we develop in this section will handle 
both the radial phase and optically-thin phase. 

The ideal picture for radial phase is very polarizing for the direct supernova flux that 
passes into the jet because of the off-radial direction is optically thin. To make the expressions 
a bit more realistic we will account for finite optical depth perpendicular to the axis in the 
following way. As defined by equation (14) in § 2, /Sgcat is the mean velocity position that 
photons incident on a jet reach for their first scattering. The perpendicular-to-axis optical 
depth from axis to jet surface at /3scat is Tos which is approximately given by equation (15) in 
§ 2. We take exp(— Tog) as the characteristic probability that an off-radially scattered photon 
escapes without scattering again. We take 1— exp(— Toff) as the probability that an off-radially 
scattered photon scatters again. We assume that photons that scatter multiple times in non- 
radial directions have on average an isotropic angular distribution and zero polarization: the 
multiple off-radial scattering should average away anisotropy and net polarization. 

Now recall that the direct supernova luminosity (emitted at wavelength A'_|_) passing into 
a jet is Ldir(A^)n(A6'). Multiplying this by 1 — exp(— r^ad) gives the luminosity that interacts 
in the jet: i.e., Ldir(A±)r2(A^)[l — exp(— Trad)]- We now assume that all this interacting 
luminosity escapes with a normalized angular distribution determined by 

c/(/Xgen) exp(-roff) + [1 - exp(-roff)] _ 
An 

i.e., a combination of Rayleigh and isotropic scattering with the relative amount of each de- 
termined by the probability of no second off-radial scattering exp(— Tofj). Thus the luminosity 
per solid angle in the direction of the Earth from a jet at wavelength A is 

Ldn{\'±MA9)[l - exp(-rrad)]:^ {^(/i) exp(-roff) + [1 - exp(-roff)]} , (38) 

47r 

where again the upper case of A'_|_ is for the near jet and the lower case for the far jet, A^ is 
determined by equation (30) (see § 3.3), and g{fi) applies to both near and far jets since this 
formula depends only on the square of the scattering angle cosine (see eq. [28] in § 3.3). If 
we multiply the luminosity per unit solid angle by dflEarth/ {D^dQEa.Tth) , where rf^Earth is a 
differential bit of solid angle at a jet subtended by the Earth and D'^dQEa.rth is the projected 
area of the Earth on the sky as seen from the jet, we get the flux at Earth from one jet: i.e., 

-^dir^Aj^ ji i(ZAC/j 

-^-j^ [1 - exp(-ri.ad)J l^(^) exp(-roff) + [1 - exp(-roff)J| 

= U{\'^)Q{A9)[1 - exp(-r,ad)] {^(/i) exp(-roff) + [1 - exp(-roff)]} , (39) 

where we have used equation (17) from § 3.1. The total flux jet flux at Earth at wavelength 
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A for the radial and optically-thin phases 

/radial/optically-thin (-^) = [fdiri^+) + /dir(^-)J fifr(A^) [1 — exp( — Tfad)] 

\g{n) exp(-roflf) + [1 - exp(-roff)] I . (40) 

Note that the 1 =F A^ factors multiplying the direct fluxes that occur in the reflection phase 
expressions (see eq. [34] and [36] in § 3.3) do not appear in equations (39) and (40): in the 
ideal radial picture all parts of the jets send light to the observer and not just the Earth-side 
parts of the inner ends. The ideal radial phase is when Toff -^ or exp(— Toff ) — > 1: recall, the 
exp(— Toff) and 1 — exp(— Tofr) factors correct for the less-than-ideal radial picture situation. 

In developing the equation (40), we are assumed the jets are not too closely aligned with 
the line of sight which would be at least approaching the occultation situation: as discussed 
in § 2, we are not concerning ourselves with the occultation situation in this paper. 

From equation (40) and again making use of our understanding of equation (28) for 
the Rayleigh phase function (see § 3.3), we find that the net Stokes Q parameter for the 
observer-scattered flux is 

<5radial/optically-thm(A) = [/dir(A+) + /dir(A'_)] rifr(A^) [1 — exp( — Trad)] 

(^) (^' - 1) eM-rosnm , (41) 

where we have again introduced the correction factor .^(A^) (see eq. [26] in § 3.3) to account 
for the range of scattering planes. Note the 1 — exp(— Tog) term in equation (40) accounts 
for unpolarized flux, and thus the x- and y-fiuxes it governs are equal and it cancels out of 

'Vradial/optically-thiri' 



3.5. General Expressions for All Phases 

The flux and Stokes Q parameter expressions of §§ 3.3 and 3.4 were developed from 
idealized limiting-case pictures. The state of the jets in the times between when these 
pictures apply is complex, but it would be useful to have general expressions that cover 
the intermediate times as well. We will construct general expressions that are valid in 
the reflection and ideal radial/optically-thin phases and act as interpolation formulae for 
the intermediate times. (The ideal radial phase is when t^s -^ recall.) We are guided 
somewhat by an intermediate phase picture. 

The reflection phase picture (see § 3.3) predicts that the fraction q of the photons 
incident on a jet undergo a single scattering and are highly polarized. The rest of the 
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photons (i.e., the fraction 1 — q) are multiply scattered and emerge unpolarized (or rather 
randomly polarized and so yielding no net polarization). The radial phase ideally assumes 
that all scattered photons are scattered as if they had undergone a single scattering event 
and is therefore a more polarizing phase ideally. Less ideally, we estimate the polarized 
fraction of the escaping scattered photons is exp(— Toff), where Tqs is the off-radial optical 
depth (i.e., the optical depth from the axis to the surface along a perpendicular to the axis) 
at the mean radial optical for the first scattering of a radial photon (i.e., at /^scat)- We will 
argue that there can be no local minimum in polarization between the reflection and radial 
phases. First, it is plausible that the ratio of singly scattered to all scattered photons should 
be roughly constant as the the reflection phase changes to the radial phase: both phases are 
optically thick in the radial direction and all photons are scattered: the photons just tend 
to penetrate more deeply as the jets expand with time. Just as the jet becomes optically 
thin perpendicular to the axis (i.e., t^s — > 0) the fraction of multiply scattered photons that 
are polarized (i.e., polarized in a non-random way) should increase from zero to all multiply 
scattered photons ideally: recall our description of the ideal radial phase scattering at the 
beginning of § 3.4. Now net polarization is set by the non-randomly polarized, scattered 
photons. Since these should increase in abundance as we go from the reflection phase to the 
ideal radial phase, we should expect polarization to increase monotonically. 

We now can derive general expressions that handle the the reflection and radial/optically- 
thin phases and interpolate between them. The discussion just above suggests replacing the 
q and 1 — g factors in the reflection phase equations for flux and Stokes Q parameter (see 
eq. [34] and [36] in § 3.3) and the exp(— Tofj) and [1 — exp(— Toff)] factors in the counterpart 
radial/optically-thin phase expressions (see eq. [40] and [41] in § 3.4) by, respectively. 



q 
1 


for Toff ^ oo; 
for Toflf -^ 
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- q for Toflf - 
for Toflf - 


* oo; 
^0, 



W = g+(l-g)exp(-Toff) = <i^ .II 1°" ,7' (42) 

and 

since the change from reflection to radial phases can plausibly be identifled as when Tofr 
drops below about 1. With these replacements the reflection phase and radial/optically-thin 
phase expressions become more alike. If we add a factor [1 — exp(— Trad)] to the reflection 
phase equations for flux and Stokes Q parameter, the expressions for the reflection and 
radial/optically-thin phases become identical, except for the 1 =F /^ factors in the reflection 
phase expressions. To complete the general expressions, we must flnd away to smoothly 
change the 1 =F /U factors into I's as the reflection phase changes to the radial/optically-thin 
phase. 
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The (1 =F Ai) factors were based on the very special, but heuristically useful, assumption 
that the inner ends of the jets could be treated as hemispheres with round sides toward the 
supernova center. We do not put much weight on this hemisphere assumption and have no 
strong guidance on how to effect the change from 1 =F /^ factors to I's as the reflection phase 
changes to the radial phase. Therefore, no elaborate procedure is warranted. One possibility 
is multiply the /x's in the 1 =F /^ factors by 

-.= 1'-^ '°^"°^^'^ (44) 

( for Toff < 1 

so that one has 1 =F w^fi instead of 1 =F ;U in the expression. The w^ factor will give smooth 
transition between the phases and at what can be regarded as the characteristic end of the 
reflection phase when Tog = 1 the w^ factor becomes a and the l^fi factors have effectively 
changed to exactly I's. 



One could, of course, devise alternative definitions of w^ that change from 1 to as the 
phases change. A very simple alternative is just to set w^ = at all times. This simple 
alternative is warranted because we do not take the hemisphere assumption too seriously 
and also because the exact value w^/i is not too important if /dir(^+) and /dir(^-) do not 
differ too much which may often be the case. Also in actual cases where polarization is high, 
fi is likely to be closer to than to 1 which would tend to make the exact value of w^ of 
little importance. 

Below we will insert w^ factors to account for the transition of 1 =F /i factors to I's, but 
we make no final specification of w^ factors. In all our demonstration calculations in § 4, we 
take 6 = 90° (i.e., /x = 0) and the value of w^ factors is irrelevant. 

Making all the changes to the expressions for the reflection and radial/optically-thin 
phases that we suggested above and replacing 1 =i= /i by 1 =F w^fi in the reflection phase 
expressions and inserting l=pw^/x in the radial/optically-thin phases appropriately, makes the 
reflection and radial/optically-thin phase expressions identical and these identical expressions 
are our general expressions for the jet-scattered flux and its Stokes Q parameter. Explicitly, 
the general expressions are 

/jct(A) = [/dir(AV)(l-«;^^) + /dir(A'_)(l + ti;^/i)] 

f]fr(A^) [1 - exp(-r,ad)] [g{f^)w + (1 - w)] (45) 



and 



Qjet(A) = [U{K){l-W,fl) + U{\'_){l+W,fl)] 

nUm [1 - exp(-rrad)] (I) (/i' - 1) w^{A9) . 



(46) 
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These expressions interpolate smoothly between the reflection phase and the ideal radial 
phase (where Tos -^ 0) and also handle the optically-thin phase, of course. This can be seen 
by comparing in detail the expressions to equations (34) and (36) in § 3.3 and equations (40) 
and (41) in § 3.4. 

For actual analysis we would like the ratio R of jet-scattered to direct flux (which we 
call the relative scattered flux) and the polarization expression for the total flux. Dividing 
equation (45) by /dir(A) gives the relative scattered flux 



R{\) 



/jet (A) ^ r /dir(AV)(l - w,fi) + /dir(A^,)(l + w,fi) 

/dir(A) [ /dir(A) 

fifr(A^) [1 - exp(-r,,d)] [g{f^)w + (1 - w)] . (47) 



For polarization we divide the absolute value of the <5jct(A) expression by total flux /dir(A)[l- 
i?(A)] and multiply by 100 % since polarization is conventionally given as a percentage: 



P(A) 



/dir(A+)(l - W^fi) + /dir(A'_)(l + W^fi) 



/dir(A)[l + i?(A)] 
f}fr(A^) [1 - exp(-rrad)] (I) (1 - /i') wC{A9) x 100 % . (48) 



Recall we assume that the direct supernova flux is unpolarized and therefore has zero Stokes 
Q parameter: thus equation (48) gives the total supernova polarization for the bipolar jet 
model. The 1/[1 + -R(A)] factor in the polarization expression is probably usually negligibly 
different from 1 for cases where our approximations are best: i.e., when jets subtend relatively 
little solid angle (i.e., when Qit.{A6) is small), and thus are relatively narrow and scatter 
relatively little flux. 

If one is trying to obtain R and P from observed flux, then one can substitute /obs(A)/[l+ 
i?(A)], /obs(A'+)/[l + i?(A'+)], and /obs(A'_)/[l + i?(A'_)] into equations (47) and (48) for U{\), 
/dir(A'_,_), and /dir(A'_), respectively. Usually one can then drop the inserted 1/(1 + i?) factors 
as being negligible and then no complex solution for i? as a function of wavelength is required. 
In fact, the inserted 1/(1 + i?) factors probably often partially cancel out of the equations: 
they exactly cancel out if R is constant with wavelength. A case where R might become large 
and where a complex solution might be needed for high accuracy is in the trough features 
of P Cygni lines where /dir(A) becomes small while [/dir(A',_) + /dir(A'_)] could stay large due 
to the redshift effect. But given that the bipolar jet model is crude, complex solutions for R 
are probably a waste of effort almost always. 

Our general expressions for R and polarization P depend on 5 independent basic param- 
eters: 6 (the angle of the jet axis from the line of sight) or /i = cos 6, A6 (the half-opening 
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angle of the jets), /3in (the inner end of the jet in velocity), b = A/9rad/An (the ratio of jet 
velocity length to inner end velocity of jet), and Trad (or if the time-delay effect is included 
tred^ see § 3.2). The r^ad parameter can be resolved into time t and the Trad.o^o parameter 
(see eq. [9] in § 2) for an analysis of the time evolution of polarization: in this case Trad.o^o 
replaces Trad as one of the 5 independent basic parameters. The q parameter, which is an 
argument of the w factor (see eq. [42] and [43] above) is also independent, but we consider its 
value as fixed at 1/4 (see § 3.3, esp. eq. [23]), and so exclude it from our list of independent 
basic parameters. The non-basic parameters in the R and P depend on the basic parameters 
as follows: A'j_ depends on 9, /3in, b, and r^-ad (see eq. [13] and [14] in § 2 and eq. [30] in § 3.3); 
w depends on A^, b, and r^ad (and on q too) (see eq. [13] and [15] in § 2 and eq. [42] and [43] 
above); w^ using our first prescription depends on A^, b, and r^ad, (see eq. [13] and [15] in 
§ 2 and eq. [44] above) and on nothing if we set it to zero. 

Given the 5 independent basic parameters, one can correct the synthetic supernova 
spectrum /dir(A) of any spherically symmetric calculation to include the additional flux from 
a jet by multiplying this spectrum by 1 -|- -R, where R itself is evaluated using the synthetic 
/dir(A). An observed spectrum (which naturally includes jet flux if there are jets) can by 
corrected to the direct supernova spectrum by dividing by 1 -|- i? or multiplying by 1 — i? 
assuming R small. In this case, R is evaluated using the observed /obs(A) (as discussed 
above) which, of course, includes jet flux if there is any: this, however, gives a direct flux 
that is first order correct in small R. 

There is considerable degeneracy among the 5 independent basic parameters and finding 
a uniquely good set for a single continuum polarization observation would be difficult. A fit 
to an observed time evolution of continuum polarization may break some of the degeneracy. 

For example, consider 6 and A^. The polarization increases with 6 up to ^'s maxi- 
mum value of 90° degrees through the 1 — /i^ factor. But the polarization also scales with 
r2fi.(A^)^(A^) which is the most obvious manifestation of the A^ parameter: f2fi.(A^)^(A^) 
increases with A^ from zero for A^ = 0° to a maximum of 0.1046 at A^ ^ 56.48° (which is 
beyond where our picture of a narrow jet has plausible validity) and then decreases with A^ 
to zero at A^ = 90°. Clearly, 6 and A^ could be varied in a compensating way while holding 
the polarization constant or nearly. Thus, there is a strong time-independent degeneracy 
between 6 and A^. 

Now A^ also affects the w factor (see eq. [42] above) through the Tos optical depth (see 
eq. [15] in § 2). Since w varies in time due to the time variation of r^ad (see eq. [13] and [15] 
in § 2), the degeneracy of A^ and 6 may be broken by fitting to an observed time evolution 
of continuum polarization. However, two other parameters come into Tos'- b (see eq. [15] in 
§ 2) and Trad.o^^o (which replaces r^-ad as an independent basic parameter in a time-dependent 
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calculation) (see eq. [9] in § 2), and these can be used to compensate variations in A^ to 
some degree. The situation is further comphcated by the fact that the time variation of Trad 
also directly affects the polarization through the 1 — exp(— Tj-ad) factor. 

Clearly, even when studying an observed time evolution of continuum polarization, 
parameter degeneracy will be a problem in finding a uniquely good set of the 5 independent 
basic parameters. Any constraints on the jet structure from non-polarimetric observations 
or theory would help to break the degeneracy. 

Another parameter degeneracy example involves the redshift effect. The redshift effect 
depends on 6 (through /i) and /3scat (see eq. [30] in § 3.3). Now /^scat depends on 6, /^m, 
and Tj-ad (see eq. [13] and [14] in § 2). The redshift effect manifests itself most prominently 
in polarization spectra: we study this manifestation in § 4.3 and 4.4. If 9, b, and Trad can 
be determined from the time evolution of the continuum polarization then that breaks the 
degeneracy with /Sm which can be determined from analysis of polarization spectra since f3in 
has a strong effect on the redshift effect (see eq. [14] in § 2 and eq. [30] in § 3.3). 

Even if a unique set of the 5 independent basic parameters cannot be obtained from 
polarimetric observations, it would still be very interesting if the bipolar jet model could be 
shown to be adequate for those observations. 



4. Demonstration Calculations 

In this section, we present demonstration calculations for continuum polarization and 
polarization spectra. In the continuum polarization case, it is of especial interest to see if the 
bipolar jet model can yield a continuum polarization time evolution that rises to a maximum 
and then declines with time (see § 4.2) and if it can yield the inverted P Cygni polarization 
profiles of lines (see §§ 4.3 and 4.4). 

For simplicity we neglect the time-delay effect (see § 3.2) in all our calculations. This 
means that we use a single reduced time tred (which specifies the radial optical depth Trad 
through eq. [10] in § 2) to calculate the polarization for each epoch or we use Trad as one 
of the independent basic parameters. Also in all our calculations, we set 6 = 90° (i.e., the 
jets are aligned perpendicular to the line of sight): recall 6 is one of the 5 independent 
basic parameters (see § 3.5). This 6 is the maximally polarizing choice as one can see from 
equation (48) in § 3.5, where the main 6 dependence is through 1 — fi"^ where n = cos6. 
(The 9 parameter also affects the polarization through the redshift effect: see equation [30] 
in 5 3.3.) 
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In § 4.1 we consider continuum polarization with varying A^ and in § 4.2 continuum 
polarization as a function of time. In § 4.3, we discuss polarization line and continuum 
features. Synthetic polarization spectrum calculations are presented in § 4.4. 



4.1. Continuum Polarization Calculations with Varying A^ 

For the strictly continuum polarization calculations in this subsection and § 4.2, we 
assume the direct flux from the bulk supernova is constant with wavelength: this is the 
neutral choice. Thus we have /dirlA'^) = /dir(A'_) = /dir(A) for these calculations and in 
equation (48) in § 3.5 the factor 

/dir(AV)(l - CJ^/i) + /dir(A'_)(l + UJ^fl) _ 

/d.(A) - • ^^'^ 

note that the values of cu^ and yU are irrelevant in the flux ratio factor for a constant continuum 
flux. The assumption of wavelength-independent direct flux effectively eliminates the redshift 
effect from the calculations (see eq. [30] in § 3.3). This means that the independent basic 
parameter /3in and the dependent parameter /5scat (see eq. [14] in § 2) never come into the 
calculations. Thus, we do not specify values of /5in in this section or in § 4.2. 

In this subsection, we are interested in seeing how continuum polarization varies with 
A^ with other parameters chosen for a high polarization case in the radial phase. To get a 
large polarization we set 6 = oo (which makes the jets infinitely long compared to distance 
from the supernova center to inner jet surface: see eq. [4] in § 2) and r^ad >> 1 (which 
means all photons entering the jet are scattered). These choices make Trad/^ = 0, Tmean = 1, 
and Toff = A6 (see eq. [13] and [15] in § 2) and 1 — exp(— Trad) ~ 1: a fairly small Tqs and 
1 — exp(— Tj-ad) nearly 1 produce high polarization (see eq. [48] and [42] in § 3.3). (As we 
will show in Figure 1 in § 4.2, slightly, but distinctly, higher continuum polarizations than 
for Trad >> 1 can be achieved for 6 = oo and A9 > 30° as Trad dechnes with time toward the 
neighborhood of 1.) 

For the specified 3 independent basic parameters (e.g., 6 = 90°, b = oo, and Trad >> 1) 
and independent basic parameter (3^^ unspecified since it is irrelevant. Table 1 shows, as a 
function of the independent basic parameter A9, the runs of Q{^{A9), C,{A9), w, R, and 
polarization P. Equation (48) in § 3.5 shows how polarization depends on the dependent 
parameters Q{,:{A6), $,{A6), w, and i? in a direct sense. The competing effects of A^ through 
the dependent parameters yield a polarization maximum of about 6.63% at A^ ^ 45.6°. 
This opening angle is beyond where our picture of a narrow jet has plausibility. However a 
continuum polarization of 4 % is reached at A^ ^ 24° which is plausibly still a sufficiently 



-33- 



narrow jet to channel non-escaping photons mostly radially as in our ideal picture of the 
radial phase discussed in § 3.4. 

The polarization result of 4% just mentioned is encouraging for the bipolar jet model 
if it turns out that polarization of this size (which is possibly the current record continuum 
polarization for a supernova: see § 1.1) is required from the model to fit observations. 
However, the result is obtained with high-polarizing choices for 9, b, and Trad- For example 
if we make b = 0, the jets become geometrically short along their axis: the radial channeling 
effect of non-escaping scattered photons disappears. The jets in this case remain in the 
reflection phase with Tqs = cxd and w = q = 1/4 (see eq. [15] in § 2 and eq. [42] in § 3.5 and 
eq. [23] in § 3.3) until they become optically thin radially. (Having w = 1/4 when the jets 
becomes radially optically thin is an unphysical limiting case for the equation [48] since most 
photons that scatter once will escape without being scattered again and therefore will not 
be depolarized by multiple scattering as w = 1/4 implies.) Much lower polarizations than 
in the b = oo case would be found since every explicit factor in equation (48) is the same as 
for the 6 = case, except that w is distinctly smaller (i.e., w = 1/4 as compared to the w 
values in Table 1) which distinctly decreases polarization and that R increases a bit which 
is because of smaller w and slightly decreases polarization. With flducial parameters used 
in Table 1, except that b is chosen to be 0, the continuum polarization reaches a maximum 
of only about 2.8 % for AO = 51°. For AO = 24°, the polarization is only 1.3 %. 

The effect of the b parameter on the time evolution of the polarization is signiflcant. 
As long as AOr^adlb (with A6 in radians, of course) is large, the off-radial escape of highly 
polarized photons, except for reflected photons the inner end of the jet, is small. Recall 
this escape is governed in the model by exp(— Xofr) (see eq. [42] and [43] in § 3.5), where 
Toff = A^(ri.ad/& + Tmcan) (s^c cq. [15] in § 2). Only when t^q becomes less than of order 1 
while Trad is still signiflcantly larger than 1 does the ideal radial phase high polarizing effect 
that was discussed in § 3.4 begin to be approached. However, if Trad is falling below 1 when 
Toff becomes of order 1, then polarization could start falling because a decreasing number 
of photons are scattered at all in the jet. Thus, we expect localized polarization maxima in 
time evolution when Tog < 1 (requiring A^r^ad/^ ^ 1), but Trad ^ 1- Our calculations bear 
this out: see § 4.2, esp. Fig. 2. 



4.2. Calculations of Continuum Polarization as a Function of Time 

In Figures 1, 2, and 3, we show the time evolution of the continuum polarization for 
representative b values and A^ values with 9 = 90°. Since we are assuming a wavelength- 
independent direct supernova flux as in § 4.1, [/dir(A'|_)(l — w^fi) + /dir(A'_)(l + u'^/i)]//dir(A) - 
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2 and /3in is unspecified since it is irrelevant for a wavelengtli-independent direct supernova 
flux. The time used is the reduced time tred = l^/^md (see § 2, esp. eq. [10]) which eliminates 
the need to specify a rrad,o^o parameter (see eq. [9] in § 2). 

Figure 1 shows the b = cases (solid lines) and b = oo cases (dashed lines) for a range 
of A^ values including 46° which yields near maximum polarization for the b = oo case (see 
Table 1). The only difference between the 6 = and b = oo cases is that in the former 
w = q = 1/4 for all times and in the latter w = q + {1 — q) exp(— A^Tmcan) (see eq. [13] 
and [15] in § 2, eq. [23] in § 3.3, and eq. [42] in § 3.5). In the b = oo cases, w will rise 
monotonically from w = q + (l — q) exp(— A6') at early times (when Trad >> 1 and Tmean = 1: 
see eq. [13] and [15] in § 2) to 1 at late times (when Tmean ^0): w in the b = oo cases will 
always be greater than q = 1/4. 

The early time behavior of the 6 = cases is the first kind of polarization plateau phase 
discussed at the end of § 2. The early time behavior of the b = oo cases is the second kind 
of polarization plateau phase likewise discussed at the end of § 2. The first kind of plateau 
phase is unavoidable and should be observable if a supernova with jets of the kind we have 
posited is observed early enough. The second kind can only exist if 6 >> 2, and is probably 
physically unrealizable for cases of significant polarization as we discuss below. 

A representative set of the Table 1 predictions for b = oo with Trad >> 1 are illustrated 
by the b = oo cases in Figure 1 at times trcd significantly less than 1. As mentioned in § 4.1, 
the b = oo cases can actually rise slightly, but distinctly, above the Table 1 predictions for 
A^ > 30° when Tj-ad declines toward the neighborhood of 1: the actual localized maxima in 
the figure occur at t^d ~ 0.54 (rrad ~ 3.5) for A^ = 30° and at trcd ~ 0.66 (rrad ~ 2.3) for 
A^ = 46°. Declining Trad tends to lower polarization by reducing the number of photons 
scattered in the jet. But, as discussed in § 4.1, the competing effect of declining Tog (which 
increases w and thus tends to increase polarization) briefly dominates and the ideal radial 
phase high polarizing effect is approached. The calculations show that this situation causes 
distinct polarization maxima to appear for sufficiently large A^. The 6 = cases are at all 
times lower in polarization than the b = oo cases of the same A^ because the 6 = case 
w factor is everywhere smaller, except all cases have zero polarization when tred = oo. The 
6 = cases decline monotonically with tj-ed'- their only time dependence comes through the 
1 — exp(— Trad) factor which declines monotonically with time (see § 3.5, eq. [48]). 

Figure 2 shows the 6=1 and 6 = 3 cases for the same A^ values as used for Figure 1. 
These cases further illustrate the effect of the 6 parameter on the time evolution. As long 
as A^rj-ad/6 >> 1 (i.e., at early times with A^ in radians, of course), the off-radial escape 
of highly polarized photons is small and the polarization is just as in the 6 = case: i.e.. 
Toff >> 1, and so w = q = 1/4. This early phase is the first kind of polarization plateau 
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phase discussed at the end of § 2. 

When A^Trad/^ becomes small the ideal radial phase effect become high and polarization 
begins to rise. However, sometime thereafter Trad becomes small and the polarization falls as 
overall scattering declines. Thus, we expect, as discussed at the end of § 4.1, a polarization 
maximum in time evolution when A^Trad/^ ^ 1 (i-e., there is high off-radial escape probability 
for polarized photons), but r^ad ^ 1 (i-e., there is still strong scattering in the jet). Figure 2 
shows this is the case. Note the polarization maxima are higher and shifted to earlier times 
for the 6 = 3 cases relative to the corresponding 6=1 cases: this is because for the 6 = 3 
cases there are longer time periods starting earlier where the conditions A^rrad/6 ^ 1 and 
Trad ^ 1 hold and a greater degree of the ideal radial phase effect is reached. Note also that 
though polarization increases with greater A^, increased A^ tends to shift the polarization 
maxima to later times because again the condition A6T^g,d/b ^ 1 has to be approximately 
met for the maxima to occur. 

The maxima in Figure 2 do not reach the high polarizations and the second kind of 
polarization plateaus obtained for corresponding 6 = oo cases (see Fig. 1) which shows that 
there is never a time when A^ri.ad/6 is negligible in the evaluation of Tofj while r^ad ^ 1- For 
instance, the 6 = 3, A^ = 30° case reaches a maximum polarization of 4.06 % at tred = 0.74 
whereas the 6 = oo, A^ = 30° case had polarization greater than 5 % for tied ^ 0.75. 

The prediction of a period of rising polarization with time given by the calculations for 
Figure 2 is qualitatively in agreement with the observations of SN 1987A (Jeffery 1991b) and 
SN 1999em (Leonard et al. 2001). The prediction of eventual decline of intrinsic polarization 
to a very low level is believed to be in agreement with observations for SN 1987A (Mendez 
1990; Jeffery 1991b). An eventual decline in supernova polarization with time must occur in 
any model depending on electron scattering polarization since the electron optical depth must 
decline eventually as the supernova ejecta expands. Our polarization expression (eq. [48] in 
§ 3.5) predicts that in all cases polarization will start to fall when Trad declines to of order 1. 

Recall ~ 4 % polarization is a possible record for continuum polarization for a core- 
collapse supernova (see § 1.1). The 6 = 3, A6 = 30° jet polarization shown in Figure 2 
reaches level. We consider bipolar jets with these parameters a physically plausible case in 
which the jets are relatively narrow. It is satisfying that a plausible bipolar jet case can 
reach the possible record polarization. 

In Figure 3, we show the polarization time evolution curves for A^ = 10° and a range 
of 6 values from to infinity. As we would expect from our discussion of the Figure 2 curves 
above, the polarization maxima in Figure 3 get higher and broader as 6 is increased since 
the time period of A^rrad/6 small while Trad ^ 1 is longer and there is a stronger ideal radial 
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phase effect. For b > 100, there is the second kind of polarization plateau. 

Physically b values over 50 are probably impossible since they give A/3rad = Pinb > 1 (i.e., 
outer end jet speeds greater than light) even for modest f3i^ = 0.02 (i.e., jet inner end speed 
of about 6000 km/s) which is of order of the velocity of supernova photospheres for much 
of the observable epoch. (Recall as A(3rad becomes more relativistic our treatment becomes 
less adequate because we have neglected relativistic effects in our formalism and because 
in the present calculations we are neglecting the time-delay effect.) Reducing A^ allows 
the second kind of polarization plateau to be reached with smaller b values. For example, 
A6 = 1° in calculations otherwise like those of Figure 3 causes the polarization plateau to 
appear for b > 10, but, as one would expect from equation (2) (§ 2) and equation (48) 
(§ 3.5), the polarization scales down everywhere by a factor of about 100 from the Figure 3 
calculations: the A^ = 1° curves have polarizations of order 0.01 % or less which are currently 
observationally insignificant. (We believe that typically the lower limit of observationally 
significant polarization is perhaps 0.1%. Depending on circumstances the lower limit could 
vary up and down from this, of course.) 

The upshot of the discussion in last paragraph is that the second kind of polarization 
plateau which is evident in Figures 1 and 3 is probably not physically realizable for high 
polarization cases. It seems unlikely that the inclusion of full relativistic effects and time- 
delay effects will ease the requirement for very high b values in order to get the second kind 
of polarization plateau for cases of significant continuum polarization. Thus the second kind 
of polarization plateau is only a mathematical prediction of the bipolar jet model for cases 
of significant continuum polarization. We do not expect it to be observed. Consequently, 
the very large and infinite b value cases seen in Figure 3 and the infinite b value cases of 
Figure 1 only illustrate mathematical trends of the bipolar jet model formalism. 

The other extreme from large b is when b is decreased toward zero (see Fig. 3 again). 
In this situation, the localized polarization maximum that can occur near trcd = 1 (i-e.. 
Trad = 1) cau first cease to be a global maximum (see the 6 = 0.1 curve) and then can 
disappear altogether (see the b = 0.01 curve) as b is reduced. 

Figure 3 also illustrates that all curves approach the same functional behavior at late 
times where they decline like Trad = ^^d (^^^ ^1- [10] i^ § 2) since 1 — exp(— r^ad) -^ Tva.d at 
late times. The late-time asymptotic curve is the same for all cases, except for the 6 = 
curve. This is because for 6 > 0, the t^s quantity goes to zero as time increases, and thus w 
goes to 1 as time increases (see eq. [13] and [15] in § 2 and eq. [42] and [48] in § 3.5). The 
6 = curve does not approach the same asymptotic curve because Tog is always infinite, and 
thus w = q = 1/4 for all times. 
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The continuum polarization calculations we have done here have been to demonstrate 
the range of polarization behaviors that might be expected from supernovae with bipolar jets 
and to illustrate the mathematical properties of the polarization formula (eq. [48] in § 3.5). 
If one had observed continuum polarization as function of time (i.e., real time, not reduced 
time), a fit to the observations could be done to determine values for 9, Ad, b, and the 
Trad^o parameter (see § 2, eq. [9]). Unfortunately, the degeneracy of the parameters which we 
discussed at the end of § 3.5 would make it difficult to determine a unique set of parameters. 



4.3. Polarization Line and Continuum Features 

Polarization spectra would have exactly constant polarization without the redshift effect 
of the jets (see § 3.3, esp. eq. [30]). With the redshift effect, P Cygni fine flux features tend 
to give rise to inverted P Cygni polarization profiles: high polarization at the blueshifted 
trough of P Cygni fiux profile and low polarization at the near-line center maximum of a 
P Cygni fiux profile. This kind of line polarization profile is also produced by other supernova 
polarization models (see § 1.1 for other supernova polarization models) and is, in fact, present 
or partially present for some observed supernova P Cygni lines (e.g., Jeffery 1991b, Wang et 
al. 2001, and Leonard et al. 2001). That the bipolar jet model can yield such a profile is a 
success. However, the observed profiles often show shifts in position angle which the strictly 
axisymmetric bipolar jet model cannot yield. Thus, the axisymmetric bipolar jet asymmetry 
even if it exists cannot be the only asymmetry in many cases and maybe in all cases. 

The inverted P Cygni profiles in the bipolar jet model arise as follows. Because of 
the redshift effect the jet-scattered, polarized flux originates in direct flux spectrum from a 
bluer part of that spectrum than the part of that spectrum to which it (i.e., the jet-scattered, 
polarized flux) is added to make up the total flux spectrum. Now if this bluer flux is redshifted 
from a continuum region to another continuum region there will be a certain scattered flux to 
direct flux ratio and a certain polarization in the second region. If the bluer flux is redshifted 
from a continuum region to P Cygni trough region, then there will be a higher scattered flux 
to direct flux ratio and higher polarization in the second region: the higher polarization is 
because there is less dilution of the scattered, polarized flux by unpolarized direct flux. If 
the bluer flux is redshifted from a continuum region to P Cygni emission peak region, then 
there will be a lower scattered flux to direct flux ratio and lower polarization in the second 
region: the lower polarization is because there is more dilution of the scattered, polarized 
flux by unpolarized direct flux. The net result of scattered continuum flux redshifted to a 
P Cygni line, is an inverted P Cygni polarization profile associated with P Cygni fiux profile: 
higher than continuum polarization in the trough region; lower than continuum polarization 
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in the emission region. 

Now P Cygni line widths are caused by Doppler shifts on scattering in the fast moving 
ejecta at and beyond the photosphere: the Doppler shifts correspond typically to velocities 
of order 1 or 2 times the photospheric velocity /9ph: the shifts can be larger for very strong 
lines. Since we assume the jets and their main scattering region at j3scat to be well outside the 
photosphere (i.e., /5scat >> /3ph), the jet redshift is usually going to be larger than the typical 
P Cygni line width. Thus, jet-scattered flux redshifted into a P Cygni line wavelength region 
will usually come from an unrelated bluer part of the spectrum which could be continuum 
or could be complicated by other P Cygni lines possibly overlapping each other. But if the 
P Cygni line is a strong one, then usually the redshifted, jet-scattered flux that is added to 
it will be more continuum- like than the P Cygni line itself. Thus for strong P Cygni lines or 
P Cygni lines sufficiently well isolated by continuum regions, one usually expects inverted 
P Cygni polarization profiles by the mechanism discussed in the last paragraph. Admittedly, 
a complex direct flux spectrum with many overlapping P Cygni will give rise to a complex 
polarization spectrum in which the inverted P Cygni profiles might not be clearly present. 

Purely continuum polarization features can arise if the jet is very fast: i.e., /^scat >> /3ph- 
This is because there will be very large redshifts and variations in the direct flux continuum 
level will affect the polarization spectrum. For example, say that the jet redshift moves 
direct flux from the ultraviolet (UV) region to the blue part of the optical and from the blue 
part to the red part of the optical. Now the UV flux of a supernova is very low compared to 
the optical in the observable epoch after very early times because of strong line-blanketing. 
In this case, the jet-scattered UV flux that is redshifted into the blue part of the optical 
will be comparatively very low and even if highly polarized might be insufficient to cause 
significant polarization of the net fiux. On the other hand, in the red part of the optical, 
the polarization should be higher than in the blue part since the jet-scattered flux from the 
blue part is shifted there. Thus, polarization should rise strongly across the optical in this 
example situation. 



4.4. Polarization Spectrum Calculations 

Using the expressions for relative flux R and polarization P (eq. [47] and [48] in § 3.5) 
and a set of the 5 independent basic parameters for the bipolar jet model, one can calculate 
the direct-plus-jet spectrum and polarization spectrum for any synthetic spectrum calculated 
in spherical symmetry. One can also take any observed flux spectrum and a set of the 5 
independent basic parameters and calculate a direct spectrum (see § 3.5 for the method) 
and a polarization spectrum. If one has an observed polarization spectrum, one can do a 
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fit to that spectrum by adjusting the parameters. Unfortunately, the effect of parameters is 
very degenerate since they combine in a complex way to control the level of polarization as 
discussed at the end of § 3.5. An analysis of the time evolution of the continuum polarization 
might allow one to at least partially break the degeneracy and obtain values for 9, A6, b, 
and Trad.o^o (s^^ § 3.5 and 4.2). An analysis of the polarization spectrum would then allow a 
determination of /9in parameter since it has a strong effect on the redshift effect (see eq. [14] 
in § 2 and eq. [30] in § 3.3). 

As a demonstration we have calculated the direct spectrum and polarization spectrum 
for an observed spectrum of SN lb 1999dn (in NGC 7714) from 17 days after V maximum 
(i.e., about 37 days after explosion assuming a rise time of about 20 days: see § 1.2). The 
observed spectrum is from Matheson et al. (2001). The photospheric velocity /3ph of the 
supernova at that epoch was 0.02 (i.e., fph = 6000 km/s) (BBK). For the 5 independent 
basic parameters we have chosen ^ = 90 % (for simplicity and high polarization), A^ = 20 % 
(for high polarization and yet still to have a relatively narrow jet: it is also the fiducial 
value from § 2), 6 = 1 (for simplicity: it is also our fiducial value from § 2), /9in = 0.06 
(for simplicity since it is exactly 3 times the photospheric velocity and also for a significant 
redshift effect), and r^ad = 0.67455 (for relatively high polarization in the radial phase of the 
jets and to obtain the fiducial mass of § 2: see below). 

The dependent parameters implied by the independent basic parameters are trcd = 
1.2176, Tmcan = 0.300, A/3rad = 0.06, Toff = 0.340, and /^scat = 0.0867. Note u^ = for 
Toff < 1 by our prescription equation (44) (§ 3.5), but its value is irrelevant since 9 = 
90% (see eq. [48] in § 3.5). If we neglect the redshift effect and set [/dir(A'_,_)(l — w^fi) + 
/dir(A'_)(l + w^/u)]//dir(A) = 2, R and polarization are, respectively, 0.0238 and 1.56%: these 
are characteristic R and polarization values for the calculation. Interpolating from Figure 2, 
we know for the specified Trad and reduced time that the continuum polarization would be 
declining from its time evolution maximum. 

If we choose /ig = 4 and t = 37 days, then equations (7) and (8) (§ 2) imply that each jet 
has mass 0.1 Mq and characteristic kinetic energy 0.724 x 10^^ ergs = 0.724 foe. This mass 
and kinetic energy were the fiducial jet values of § 2 and our choice of precisely 0.67455 for 
Trad was made to yield them. 

Figure 4 shows the observed flux spectrum and the calculated direct flux spectrum. 
The direct spectrum was calculated using a R value which in turn was calculated from the 
observed fiux by the method suggested in § 3.5. Since the characteristic R value for the 
calculation is small (0.0238: see above), this method is quite accurate. The scale of the 
spectra is arbitrary. The spectra have been corrected for the redshift due to the host galaxy 
heliocentric velocity (2798 km/s) and for the host galaxy Galactic foreground reddening 
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(i?B-v=0.052) using the reddening law of Cardelli et al. (1989) with the revision of O'Donnell 
(1994). Both host galaxy correction parameters were obtained from the NED database 
(http : //nedwww . ipac . caltech . edu/). 

Since the direct spectrum is on average only about 2.4% lower than the observed spec- 
trum, the two spectra virtually overlap. The effect of the jet-scattered flux on the observed 
spectrum is probably negligible compared to observational error. Having a negligible effect 
on an observed spectrum at the same time as producing significant polarization was, of 
course, a main desideratum from the bipolar jet model since this desideratum would resolve 
the potential paradox discussed in § 1.2. 

Figure 5 shows the calculated bipolar jet model polarization spectrum overlaid on the 
observed flux spectrum. Note that the polarization spectrum can only extend as far to 
the blue as the jet-redshifted blue end of the flux spectrum. Inverted P Cygni polarization 
profiles associated with flux spectrum P Cygni lines are evident: the one associated with 
the H and K Ca II P Cygni line (trough centered at ~ 3750 A) is the most clear example. 
The line blending in the flux spectrum, of course, somewhat distorts the inverted P Cygni 
polarization profiles. The largest polarization maxima, not counting the far red end of the 
spectrum, reach over 4 %. The complex high polarization features at the red end of the 
spectrum rise higher: the highest is about 9 %. We cut off these high polarization features 
in the figure to preserve the clarity of the rest of the polarization spectrum. In any case, the 
red end of the flux spectrum from which the these high polarization features are partially 
determined is probably uncertain. 

For Figure 6, we have changed the independent basic parameter /3in to 0.12. The changed 
dependent parameters are A/9rad = 0.12 and /Sgcat = 0.173. The implied mass of each jet 
is changed from the 0.1 M© to OAMq: this follows from holding the ratio M/f3f^ fixed in 
equation (7) (see § 2) so that the other values in that equation can stay unchanged. This 
mass is a relatively large fraction of a typical SN lb where total ejecta mass is perhaps only 
4Mq or less (e.g., Hachisu et al. 1991). 

In the Figure 6 calculation, the wavelength of jet-scattered flux is increased by a factor 
of 1.173 (see § 3.3, eq. [29] with yU. = 0) which is so large that the jet-scattered flux comes from 
direct flux spectrum wavelength regions where the continuum levels can be rather different 
that those of the wavelength regions to which the scattered flux is added to give the net 
spectrum. A result of this large redshift is that polarization at the blue end of the optical 
is somewhat reduced in overall continuum trend relative the Figure 5 case. The polarized 
flux in the blue part originates in even bluer flux of the direct spectrum: but since the direct 
spectrum declines blueward on average, there is less polarized flux at a given wavelength 
than if the redshift were smaller. 
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Figure 6 like Figure 5 shows inverted P Cygni polarization profiles. 



5. Conclusions and Final Discussion 

We have proposed a bipolar jet model for supernova polarization. The motivation for 
this model is to resolve the potential paradox set by BBK analysis which suggests strong 
spherical symmetry for Type lb supernovae (see § 1.2) and the expectation that Type lb 
supernovae should be highly polarized in the continuum up to perhaps ~ 4 % (see § 1.1). For 
the model, we have derived simple, non-relativistic approximate analytic formulae for relative 
scattered flux R (i.e., for the ratio of scattered flux to direct flux from the bulk supernova) 
and for polarization (see eq. [47] and [48] of § 3.5). The polarization formula qualitatively 
reproduces two main features observed for at least some supernovae: (1) the rise of continuum 
polarization to a maximum and then a decline with time (see § 4.2); (2) the inverted P Cygni 
polarization profiles of lines (see § 4.4). The level of continuum and line polarization obtained 
reaches that observed for supernovae (i.e., of order a percent or a few percent: see § 1.1) for 
plausible jet model parameters: in fact, one set of plausible parameters does allow continuum 
polarization to reach ~ 4% (see § 4.2). The position of angle polarization is a constant for 
the bipolar jet model since it is a strictly axisymmetric model and the model does not allow 
90° shifts in position angle which are possible for axisymmetric models in general. Since 
position angle variation is observed, the model cannot account for all the polarization of all 
supernovae. 

The bipolar jet model predicts that two kinds of the polarization plateau phases in time 
can exist (see § 2). The first kind of plateau phase is unavoidable and should be observable if 
a supernova with jets (of the kind we have posited) is observed early enough: we show cases 
of this plateau phase in Figures 1, 2, and 3 in § 4.2. The second kind can only exist if 6 >> 2 
(see § 2): we show cases with and without the second kind of plateau phase in Figures 1, 2 
and 3 in § 4.2. The second kind is probably physically unrealizable for jets wide enough to 
give observationally significant continuum polarization (see § 4.2). We have not yet examined 
existing supernova polarization data to find evidence for the kinds of polarization plateaus 
predicted by the bipolar jet model. 

We derived the analytic formulae for relative scattered fiux and polarization using very 
simplifying assumptions. We have treated the electron optical depth structure of the jets very 
crudely and neglected any line or other opacity in them. The geometry of the supernova-jet 
structure has also been treated crudely (see § 2). The time-delay effect (see § 3.2) can be 
included in the model, but it can be neglected for simplicity as we have done in deriving 
the analytic formulae and in the demonstration calculations. We have assumed the bulk 
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supernova (i.e., the supernova not counting the jets) is spherically symmetric: this cannot 
be completely true although to resolve the potential paradox it must be roughly true. 

There are 5 independent basic parameters for the bipolar jet model: 9 (or ^ = cos 9), 
^9, Trad, All! and h (see § 3.5) which unfortunately are rather degenerate particularly for 
the analysis of a single polarization epoch. For an analysis of the time evolution of the 
polarization the Trad^o parameter (see § 2, esp. eq. [9]) can replace r^d as one of the 5 
independent basic parameters. In an analysis in which the time-delay effect is included, trcd 
(see § 3.2) replaces Trad as one of the 5 independent basic parameters. 

We have neglected relativistic effects in our formalism. If Acat becomes sufficiently large 
(e.g., > 0.1), then this neglect could be a significant limitation of the formalism since rel- 
ativistic effects will be important. These effects could significantly reduce the calculated 
supernova polarization relative cases that neglect them, but have otherwise the same pa- 
rameters. The essential reason is that energy is lost from scattered flux by Doppler shifting 
into the jet frame that is not regained Doppler shifting back to the rest frame of the super- 
nova center. If the jet-scattered flux decreases because of this process, then the amount of 
polarized flux decreases and net polarization decreases. Thus, our formulation may overes- 
timate polarization for jets that are too relativistic. In future work we intend to investigate 
relativistic effects in the bipolar jet model. 

Because of our simplifying assumptions, the derived formulae can only be expected to 
semi-quantitatively reproduce the effects of bipolar jets and, of course, the parameters used 
in the formulae will have only a crude relation to the parameters of actual bipolar jets which 
we expect to have rather complex structure. The formulae should be tested by Monte Carlo 
calculations. Such tests would probably show how to adjust the formulae to improve their 
quantitative accuracy. In any case, variations on the formulae we have presented are clearly 
possible. In the analysis of actual supernova polarimetric data, some variations may turn 
out to be obviously more appropriate than formulae we have given. 

The advantage of the analytic formulae is that they allow a straightforward, intelligible 
analysis of polarimetry data by fltting the data and they allow a simple flrst-order approach 
to testing the adequacy of the bipolar jet model. We have verifled the adequacy of the model 
to some degree in this paper: see the flrst paragraph of this section. If the model's adequacy 
is further verifled, then more a more elaborate testing procedure for the model should be 
considered. For example the calculation of polarization spectra using Monte Carlo radiative 
transfer for a bipolar jet model whose structure is produced — if it is feasible — by a realistic 
hydrodynamic calculation. 

In this paper, we have only considered the bipolar jet model with SNe lb in mind. The 
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bipolar jet model may apply generally to core-collapse supernovae. On the other hand it may 
only to apply to some: i.e., some SNe lb and, perhaps, some of their smaller-envelope kin. 
Type Ic supernovae. In future work we hope to test the bipolar jet model using polarimetry 
and spectropolarimetry from all types of core-collapse supernovae. 

As final word, one must emphasize that neither the bipolar jet model nor the other 
models of supernova asymmetry (see § 1.1) are exclusive. Several kinds of asymmetry can 
be present for a single supernova. For example. Type Ic supernova SN 2002ap shows a 
component of polarization that can modeled by a jet (Kawabata et al. 2002; Leonard et 
al. 2002). There is certainly another polarization component present though in SN 2002ap 
perhaps due to a inner ejecta asymmetry of order 10 % in axis ratio projected on the sky 
(Kawabata et al. 2002; Wang et al. 2003). The late direct HST imagining of SN 1987A 
confirms that in one supernova at least, there is a large inner ejecta asymmetry: the major 
to minor axis radio of images is of order 1.5/1 (Wang et al. 2002). It may be that such 
inner ejecta asymmetry (caused perhaps by non-relativistic, never-emergent bipolar jets: 
Hoflich et al. 2001) is the main kind of SN 1987A asymmetry. However, a re-analysis of early 
SN 1987A speckle interferometry (Nisenson & Papaliolios 1999) suggests that relativistic jets 
(with perhaps {3 ~ 0.8) did emerge from SN 1987A in addition to other asymmetries. These 
suggested jets, however, may have been too relativistic to have contributed significantly to 
SN 1987A's polarization: this a subject for further study. 

I thank David Branch and Eddie Baron for stimulating this work and their comments on 
it. Support for the paper was provided by the Departments of Physics of New Mexico Tech 
and the University of Nevada, Las Vegas, and the Department of Physics and Astronomy of 
the University of Oklahoma. This research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute 
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Fig. 1. — Continuum polarization as a function of reduced time for a range of A^ values 
as specified on the figure, and 6 = and b = oo. We have set 6 = 90° and the direct 
supernova flux is wavelength independent. The later condition means that [/dir(A'|_)(l — 
'^M/^)+/dir(A'_)(l+w^/u)]//dir(A)] = 2, and the /3in parameter is not used and is left unspecified. 
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Fig. 2. — The same as Figure 1, except the polarization curves are for 6=1 and 6 = 3. 
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Fig. 3. — The same as Figure 1, except that A^ is set to 10° and a range of b values as 
specified on the figure are used. 
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Fig. 4. — The observed flux spectrum and bipolar jet model direct flux spectrum of 
SN 1999dn for day 17 past V maximum. The scale of the spectra is arbitrary. The bipolar 
jet model parameters are specified in the text (see § 4.4). Both spectra have been corrected 
for heliocentric velocity redshift and foreground Galactic extinction. 
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Fig. 5. — The bipolar jet model polarization spectrum for the observed flux spectrum of 
Fig. 4 which is shown overlaid. The bipolar jet model parameters are the same as for Fig. 4 
and are specified in the text (see § 4.4). 
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Fig. 6. — The bipolar jet model polarization spectrum for the observed flux spectrum of 
Fig. 4 which is shown overlaid. The bipolar jet model parameters are as for Figs. 4 and 5 
(see text § 4.4), except the /3in parameter has been doubled from 0.06 to 0.12. This greatly 
increases the redshift effect of the jets. 
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Table 1. Runs of High Continuum Polarization and Ancillary Quantities as Functions A^ 

for the Radial Phase of the Bipolar Jet Model 



Ae f7f,,(A6») g(A6») w R P 

n (%) 











1 


1 








1 


7.6 


X 10-5 


0.9998 


0.9870 


1.1 X 10-"* 


0.0113 


3 


6.9 


X 10-* 


0.9982 


0.9617 


0.0010 


0.0986 


7 




0.0037 


0.9901 


0.9137 


0.0058 


0.5029 


10 




0.0076 


0.9798 


0.8799 


0.0119 


0.9708 


15 




0.0170 


0.9549 


0.8272 


0.0270 


1.9657 


20 




0.0302 


0.9207 


0.7790 


0.0486 


3.0939 


24 




0.0432 


0.8871 


0.7433 


0.0704 


3.9944 


30 




0.0670 


0.8270 


0.6943 


0.1107 


5.1942 


40 




0.1170 


0.7053 


0.6231 


0.1975 


6.4399 


46 




0.1527 


0.6224 


0.5860 


0.2606 


6.6262'' 


50 




0.1786 


0.5643 


0.5634 


0.3069 


6.5166 


60 




0.2500 


0.4135 


0.5132 


0.4359 


5.5421 



"Maximum continuum polarization (for the chosen parame- 
ters) as a function of A9. To better numerical accuracy it is 

6.6271108% at 45.600748°. 

Note. — The predictions were calculated using equation (47) 
for R and equation (48) for P assuming a wavelength indepen- 
dent direct flux from the bulk supernova. The fixed parameter 
values for the predictions are 6 = 90°, b = oo, and T^ad >> 1 
(see text, § 4.1). The /?;„ parameter is not specified since it is 
irrelevant to the calculation. 



